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0.  Abstract 

We  present  a  model  that  provides  predictions  for  the  occurrence  of  disruptive  deviations 
of  signal  line  activity  in — and  consequent  malfunction  of— a  microcontroller  (pC) 
subjected  to  external  irradiation  by  a  narrowband  electromagnetic  (EM)  pulse.  In  our 
model,  the  state  of  a  pC  is  completely  specified  by  giving,  for  each  of  its  relevant  signal 
lines,  the  signal  pulse  train  (SPT)  time  history  on  the  line  during  any  fixed  but  arbitrary 
time  window  of  interest.  The  occurrence  of  such  disruptive  deviations  is  observed 
experimentally  to  behave  stochastically  in  at  least  some  EM  pulse  frequency  regimes — 
for  example,  in  the  radio  frequency  (RF)  regime — and  our  model  provides  predictions  for 
the  probability  of  such  disruptions  based  upon  all  relevant  characteristics  of  the  EM  pulse 
and  of  the  pC  SPT’s.  In  the  present  paper  we  focus  our  attention  on  signals  traversing  a 
single  pC  signal  line. 

I.  Introduction 

In  this  paper,  we  are  interested  in  the  following  situation.  Initially,  an  electromagnetic 
(EM)  pulse  impinges  upon  a  microcontroller  (pC)  in  which  at  least  one  signal  line  is  active, 
this  activity  resulting  from  (1)  a  single  software  application  executing  on  the  |uC  or  from 
(2)  other  dynamic  actions  of  the  piC — “system  tasks” — e.g.,  system  clock  pulse 
generation,  program  counter  increment,  instruction  fetch,  or  the  hardware-implemented 
portion  of  interrupt  handling.  This  pulse  electromagnetically  couples  to  the  jliC  at  one  or 
more  of  its  spatial  locations  (e.g.,  at  |uC  control  or  data  lines  or  at  register  or  memory 
structures  or  at  power  distribution  system  components).  This  coupling  quantitatively 
depends  upon  the  hardware  details  of  the  jliC  and  may  cause  sufficient  deviations  of  the 
normal  time-dependent  signal  pulse  trains  (SPT’s)  on  some  of  the  signal  lines,  during 
some  fixed  but  arbitrary  finite  time  window  [0,  TE)  of  interest,  to  disrupt  proper  |uC 
functionality  (i.e.,  proper  software  application  task  or  system  task  execution)  sometime 
during  that  window.  We  term  this  disruption  jliC  “upset".  This  upset  is  observed 

experimentally  to  occur  stochastically  in  the  radio  frequency  (RF)  regime — and  we 
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present  here  a  signal-centric  model  that  provides  predictions  of  the  probability,  pupse[  of 
such  upset  based  upon  all  relevant  characteristics  of  the  EM  pulse  and  of  the  jliC  SPT’s. 

There  has  been  significant  research  into  the  effects  of  RF  signals  at  the  device  level 
(Richardson  1 979,  Firestone  2004,  Estep  201 0,  201 1 ,  Estep  et  al  201 1 ),  as  well  as  at  the 
circuit  and  chip  level  (Laurin  et  al  1 995,  Liu  &  Ho  1 995).  More  recently  the  2001  -2006  RF 
Effects  MURI  (Granatstein  2007)  provided  significant  advances  in  understanding  and 
predicting  RF  effects  (Firestone  et  al  2006,  Hemmady  et  al  2005,  Yang  &  Kollman  2006, 
Wang  et  al  2006).  Building  on  this  work  to  develop  a  system-level  model  for  RF  effects 
on  a  full  digital  system  such  as  a  PC  is  daunting,  in  part  because  of  the  vast  number  of 
transistors  that  make  up  a  modern  CPU  or  memory  module,  but  also  because  a  PC 
consists  of  multiple  printed  circuit  boards  each  containing  many  integrated  circuits, 
discrete  components  and  interconnects.  These  boards  are  connected  together,  housed 
in  a  metal  box  and  connected  to  various  peripheral  devices.  For  these  reasons  a 
microcontroller,  being  a  simple  but  complete  computer  on  a  single  chip  and  thus 
representing  an  intermediate  level  of  complexity  between  an  individual  CMOS  device  and 
a  complete  digital  system,  is  an  ideal  device  to  use  as  a  test-bed  in  attempting  to 
understand  and  predict  RF  effects  on  digital  electronic  systems  (Taylor  201 1 ,  Clarke  et  al 
2011,  Clarke  et  al  2012,  Henderson  et  al  2012).  Earlier  research  on  quantifying  the 
immunity  of  microcontrollers  to  ESD  pulses  (voltage  spikes)  (Vick  &  Habiger  1997, 
Wendsche  et  al  1 999)  indicated  that  the  susceptibility  varied  with  the  particular  instruction 
(or  micro-instruction)  being  executed  when  the  ESD  pulse  was  applied.  This  suggested 
that  carefully  timed  pulses  relative  to  clock  signals  could  yield  useful  information  about 
the  underlying  mechanisms  of  RF  interference,  and  provided  much  of  the  motivation  for 
the  current  work. 

The  remainder  of  this  paper  proceeds  as  follows.  In  Section  II  we  present  the  conceptual 
structure  upon  which  our  model  rests,  discussing  the  elements  of  a  completely  specified 
bridge  between  the  application/system  tasks  executing  on  the  |uC  and  the  resulting  normal 
SPT’s  (i.e.,  those  in  the  absence  of  an  externally  impingent  EM  pulse)  on  the  |uC  lines 
corresponding  to  the  execution  of  those  tasks.  We  also  discuss  non-normal  SPT’s  (i.e., 
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those  in  the  presence  of  an  externally  impingent  EM  pulse)  on  the  piC  lines  resulting  from 
the  normal  SPT’s  via  exposure  of  the  |uC  to  an  EM  pulse.  In  addition,  we  discuss  the 
approximations  that  we  make  in  our  treatment.  In  Section  III,  we  consider  the  interaction 
between  an  idealized,  normal  SPT,  n71,  being  carried  on  a  single  signal  line,  and  a  single 

EM  pulse,  carefully  characterizing  for  our  purposes  both  of  these  entities  as  well  as 
their  interaction  itself.  This  leads  in  Section  IV,  to  an  expression  for  the  probability  of  a 
^■EM-induced  disruptive  deviation  of  n71.  Sections  III  and  IV  are  concerned  with  only  a 

single  signal  line.  Those  results  may  be  extended  to  an  arbitrary  finite  number,  L,  of  signal 
lines  carrying  an  /.-tuple  of  SPT’s  interacting  with  one  and  the  same  EM  pulse,  but  we 
defer  our  consideration  of  multi-signal-lines  to  a  subsequent  paper.  Finally,  since  in 
practice  it  will  be  rarely  known  with  certainty  which  particular  SPT  of  those  possible  is  on 
the  line  during  [0,  TE),  we  consider,  in  Section  V,  random  (stochastic)  mixes  of  possible 
SPT’s  on  a  single  signal  line.  In  Section  VI  we  conclude. 

II.  Conceptual  Overview 

In  this  section  we  discuss  the  conceptual  structure  of  our  model  as  well  as  the  modeling 
approximations  that  we  make.  The  goal  here  is  to  build  a  conceptual  bridge  between  the 
executing  application/system  tasks  and  the  SPT’s  resulting  therefrom,  both  in  the  normal 
and  non-normal  cases.  We  understand  that  some  of  the  theoretical  entities  we  construct 
may  not  be  presently  readily  obtainable  but  there  is  nevertheless  some  merit  in  laying 
down  the  logical  structure  of  our  model. 

11.1  pC  Tasks 

As  indicated  in  the  Introduction,  we  are  interested  in  the  upset  of  a  jliC  when  a  it  is 
executing  a  ( single )  software  application  concurrently  with  system  tasks  during  some  time 
interval  [0,  TE)  of  interest.  For  us,  an  application  is  a  finite  set  of  assembly  language  (AL) 
instructions,  together  with  the  specification  of  the  order  of  execution  of — i.e.,  of  the  time 
sequence  of — those  AL  instructions.  In  general,  a  given  application  in  fact  comprises 
many  possible  AL  instruction  execution  sequences,  which  sequences  we  term 
realizations,  as  a  consequence  of  the  possibility  of  branching.  An  application  having  no 
possibility  of  branching  (including  no  possibility  of  interrupts)  has  only  one  realization.  We 
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consider  an  interrupt  request  signal  and  the  subsequent  hardware-implemented  interrupt 
handling  response  thereto  to  be  system  tasks  but,  in  contrast,  we  consider  the  execution 
of  the  stored-in-memory  interrupt  servicing  AL  instructions  themselves  to  be  part  of 
application  execution.  More  generally,  we  use  the  term  “system  task”  to  refer  to  a  |uC 
action  that  requires  no  AL  instructions  in  order  to  execute  (i.e.,  is  implemented  in 
hardware)  and  we  use  the  term  “application  task”  to  refer  to  a  pC  action  that  does  require 
AL  instructions  in  order  to  execute. 

To  facilitate  further  discussion,  we  introduce  some  constructs  that  we  will  subsequently 
employ  to  describe  |uC  tasks.  To  begin,  we  label  each  software  application  that  may 

execute  on  the  jliC  by  its  own  “application  label”,  generically  a,  and  denote  by  ?4  the  finite 
set  of  all  application  labels  arising  from  all  those  jliC  applications,  so  that  #  e  rf.  (The  null 
application,  i.e.,  no  application  actively  executing,  is  denoted  by  the  unique  label  a  =  0.) 

We  denote  by  k*\,  k%  ...,  *V,  N*>  1  if  <z  *  0,  the  finite  number  of  all  possible  realizations 
of  application  where  each  k*  is  a  finite  sequence — 

k*j  =  <i*  -| ,  i*2,  ■■■,  1  if  0 — in  the  set  $  =  {^,  ...,  iNa),  Ntj  >  1 ,  of  all  assembly 

language  instructions  available  to  the  jliC;  and  we  write  (when  a*  0)  <E^  = 

{k* ,  k%  ...,  AVI  ^  u;_mr[°’  T^N  ^ancl  we  set  ^0  =  0)-  Here  ^raxt0’  Te)  is  the  maximum 

ALI 

number  of  AL  instructions  executable  in  [0,  TE),  namely  A/TOX[0,  TE)  <  [7e/tCji/2]  +  1 ,  where 
tc  >  0  is  the  system  clock  period  with  xc  1/2  =  xc/2  and  where  [•]  denotes  the  greatest 
integer  function,  with  xc1/2  being  the  lower  limit  of  our  time  scale  resolution.  Now  each 
instruction  i*(  executes  on  the  piC  during  some  time  interval  (t *e))  for  which  (x* 

e)f  -  (x*  ()0  is,  under  normal  execution,  an  integral  multiple  of  xc1/2.  Because  some  AL 
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instructions  require  more  clock  periods  to  fully  execute  than  do  other  AL  instructions,  it  is 
necessary  for  our  purposes  to  append,  to  each  AL  instruction  in  the  ^sequence,  the  initial 

and  final  time  instants  of  the  time  interval  in  which  the  instruction  executes.  We  thus 
specify  a  finite  sequence  of  instructions  executing  during  [0,  TE)  as  the  finite  sequence  of 

ordered  pairs  [A*]#  =  <  <i«  r  [(x«  ,)0,  (x*  e)f>  >  ^  where  e  *  and  t(x</,  ,)0>  (xy,  X)  £  [°> 
Te)\  and  we  write  as  well  =  {[^]#,  [ha2 ]#,  ■  ■■,  LVf}  ^ 

A/AL I  a  1  # 

Uw:r  [4x[0,  Te)]  (and  <£0  =  0).  Note  that  (x«,),-  (x«,)0  =  n^xc,1/2  where  e 

ALI-CYC  ALI-CYC 

{1,2,...,  A/tox  },  with  A/tox  being  the  maximum  number  of  clock  V2-cycles  required 

ALI-CYC 

by  any  AL  instruction  (e.g.,  A/rrax  =  1,  2,  4,  8,  16,  24,  depending  upon  the  particular 
(.iC). 

Next,  we  write  the  finite  set  of  all  system  tasks  as  0  =  {<90,  0V  ...,  0N&),  N&  >  1  (with  0O 
denoting  the  null  system  task,  i.e.,  no  system  task  executing),  and  any  sequence  of 
system  tasks  from  0  as  dk  =  <6k  A,  0k  2,  ■■■,  0k  N*k >,  tXk >  1 ;  further,  we  write 

/§>  =  {dv  dNJ,  N&  >  1,  for  the  finite  set  of  all  such  sequences,  with  k  c  Te>  &N, 

where  A/^JO,  TE)  is  the  maximum  number  of  system  tasks  executable  in  [0,  TE).  (Note 
that  it  is  not  in  general  true  that  A/^TX[0,  TE)  <  ITe/xc  V2]  +  1  since  several  system  tasks 
may  execute  simultaneously  during  any  given  clock  Vfe-cycle.)  As  in  the  case  of  AL 
instructions,  each  system  task  0k  m  executes  on  the  jliC  during  some  time  interval 

[(xfc  rt) o>  (Tfc  n)t>  f°r  which  (x^  Jf-  (xJfc  m)0  is,  under  normal  execution,  an  integral  multiple 

of  xc>1/2  (or  is  taken  by  us  to  be  so  since  xc  1/2  is  the  lower  limit  of  our  time  scale  resolution); 
further,  since  interrupts  in  general  occur  irregularly  in  time,  it  is  necessary  for  our  purposes 
to  append,  to  each  system  task  0k  m  in  the  dk  sequence,  the  initial  and  final  time  instants 

of  the  time  interval  in  which  it  executes.  We  thus  specify  a  finite  sequence  of  system  tasks 
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executing  during  [0,  TE)  as  the  finite  sequence  of  ordered  pairs  [<*J#  =  <  <0k  m,  [(xdk  J0,  (t* 
where  0k,me@  and  [(xl  J0,  (xdk  J)  <= 

M  mST  [0  T \  ly  i 

[0,  7E);  and  we  write  as  well  4  =  {[^]#,  [<y#,  [^]#}  c  l)N™*x  ’  E  [®x  [0,  rE)]  Note 

ST-CYC  ST-CYC 

that  (4  J,  -  (^  Jo  =  rrk  mxc>1/2,  where  {1,2,...,  A/TOX  }  with  Nmx  being  the  maximum 
number  of  clock  V2-cycles  required  by  any  individual  system  task  (typically 
'7U=1)-  Further,  in  addition  to  possible  overlap  in  time  of  system  tasks,  i.e., 

[K  Jo-  K  m ))  n  [K  Jo-  (xfc  m')f)  ^  0  for  some  m,  m'  €  {1 ,  ...,  Yj  with  m  *  m',  it  is  also 

possible  that  an  AL  instruction  and  a  system  task  may  execute  simultaneously  at  the  V2- 
cycle  time  scale  (e.g.,  in  the  case  of  a  pre-fetch  capability  in  some  pC's)  so  that  [(xyj0, 

n  IK  Jo-  K  J ))  *  0  for  some  1=  1 ,  ...,  N]  and  m  =  1 ,  ...,  N\. 


We  observe  that  the  system  task  sequence  fundamentally  differs  from  the  AL  instruction 
sequence  in  that  the  latter  is  “quasi-deterministic”  in  time  while  the  former — if  it  contains 
interrupts — may  be  irregular  in  time.  That  is,  in  the  absence  of  interrupts,  the  AL 
instruction  execution  sequence  timing  is  predictable  for  any  given  application  realization 
while,  on  the  other  hand — for  either  the  active  or  quiescent  state  of  any  executing 
realization  of  some  application — the  interrupt  sequence  timing  is  in  general  not 
predictable  since,  for  example,  external  interrupts  are  activated  by  conditions  appearing 
irregularly  on  devices  external  to  the  pC  (which  devices  are  monitored  by  the  pC’s  input 
pins).  This  renders  as  non-deterministically  predictable  the  particular  realization,  h*,  of  an 

application  a,  that  is  actually  executed  in  [0,  TE)  in  any  particular  instance  of  the  execution 

of  <1,  since  ha  itself  comprises  the  interrupt-servicing  AL  instructions;  hence  the  inference 

from  any  particular  k*  instruction  stream  of  the  induced  SPT’s  on  the  |uC  signal  lines  as  a 
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result  of  the  execution  of  a  also  suffers  from  this  unpredictability.  For  us,  the  statement 


that  application  a  is  executing  during  [0,  TE)  is  tantamount  to  the  statement  that  a 
particular  realization  of  a,  say  \y  e  is  executing  during  [0,  TE).  However,  to  analyze 

the  jliO  SPT  behavior  during  execution  of  we  must  resort  to  analyzing  some  statistical 

mix  of  h*' s  that  are  possibly  executing  during  that  [0,  TE),  hence  some  statistical  mix  of 

SPT’s  that  are  possibly  executing  during  that  [0,  TE).  We  address  this  state-of-affairs  in 
Section  V  where  we  discuss  random  (i.e.,  stochastic)  collections  of  SPT’s. 

As  per  the  above  considerations,  we  may  indicate  any  possible  application-^/system  task 

combination  executing  on  the  pC  during  [0,  TE)  by  the  pair  <[k*f,  [d]#  >  e  <F£aX/§>#. 

There  is  however  an  additional  factor  that  must  be  taken  into  account:  the  time-dependent 
data  (numbers)  occupying  the  various  pC  registers  and  memory  locations.  These  values, 
being  represented  by  binary  0’s  and  1’s,  determine  which  data  signal  lines  are  low  and 
which  are  high  at  any  particular  instant,  i.e.,  the  SPT’s  on  these  data  lines.  These  values 
arise  from  user  inputs  to  the  application,  or  from  inputs  to  the  application  by  the 
aforementioned  external  devices  connected  to  the  pC,  or  as  values  computed  via 
arithmetic  manipulations  by  the  application.  In  principle  these  data  may  be  represented 
by  a  time  sequence  of  the  data  contained  in  each  register  and  memory  location  used  by 
the  application.  If  we  label  the  registers  used  by  realization  ¥) of 

application  <z  as  9i  =  {pv  p2,  ...,  pN{^},  NJftp  >1  if  <z  *  0,  and  if  we  label  the  memory 
locations  used  by  application  4- as  M  =  {pv  p2,  ...,  pN  (^},  NJX)  -  1  if  0,  then  we  may 

N  (h^)  N  (h^)  [0,  T) 

represent  the  data  history  as  <  <dp  (nxCiV2)>p^\J,  <^(^,1/2 )>q=  \J  >n  _o  E  where 

dp  {mc> i/2)  is  the  value  of  the  data  in  register  pp  at  time  mc1/2  (the  beginning  of  !/2-cycle  n), 
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CYC 

and  similarly  for  (mc  1/2),  and  N  [0,  TE )  =  [7e/tc  1/2].  Alternatively,  we  may  represent 

the  data  histories  of  the  registers  and  memory  locations  in  a  fashion  that  mimics  the 
representation  method  that  we  used  for  [haf  and  [<s]#.  To  this  end  we  note  that  since,  for 

any  given  register,  pp,  there  are  only  a  finite  number  of  instants  at  which  the  data  value  in 

that  register  actually  changes,  say  at  instants  0  <  x\p  <  ...  <  xNp  <  TE  (with  xP0p  =  0  and 

xNp  +1=  Te)  then  we  may  write  the  time  history  of  the  data  in  register  pp  as  the  sequence 


<  <cIpKp- i)>  l\P-v  \P)>  >^p+1  where  dp(xPp)  *  dp{xp_,)  and  if 


t,  t'  e  [iPup_v  xPp)  then  dp  (f)  =  dp  (t’).  We  then  write 


=  «<dPp{  for  the  histories  °f  al1  the  registers  employed 

by  Analogously,  for  the  memory  locations  employed  by  ¥*  we  write  cTm^  = 

«<dp{xq^),  [xPq_v  \q)>\Pq^>qL‘{'^\  and  we  then  write  M(k«)  = 

Finally,  we  denote  the  set  of  all  possible  data  histories  of  the  form  4(,  as  D^so  that 


d- <31,  M$j)  £  T)  a  ■ 


We  may  now  indicate  any  possible  application-^/system  task  combination  executing  on 
the  jliC  during  [0,  TE)  by  the  triple  %  =  <[h*f,  [d]#,  cf#}l m{^)>  e  ^x/xD*,,,  which  triple 


specifies  collectively  an  AL  instruction  stream,  a  system  task  stream,  and  streams  of  all 
applicable  register  and  memory  location  contents,  with  all  of  these  streams  having  timing 

information  precisely  specified.  Since  the  three  components  of  %  are  each  precisely 

specified  then  the  SPT’s  resulting  therefrom  on  all  jliC  signal  and  data  lines  are  then  in 
principle  also  precisely  inferable.  This  conversion  from  %  to  the  resulting  induced  SPT’s 
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is  the  bridge  to  which  we  have  alluded  previously.  However,  since  having  such  precise 
knowledge  is  not  realistic — because,  for  example,  of  the  randomness  in  the  timing  of  the 

interrupts  of  [<if  as  well  as  the  randomness  in  cC^M{ha)  of  the  values  and  timing  of  the 

data  sent  to  the  |uC  by  external  devices  connected  to  it.  Nevertheless,  in  order  to  proceed 
we  will  assume  initially  that  a  (precisely  specified  as  above)  %  is  given  so  that  its 

associated  resulting  SPT’s  on  all  signal  lines  are  known  precisely  as  well  (although  we 
will  not  provide  a  prescription  for  the  explicit  SPT’s  so  resulting).  For  our  purposes  it  is 

sufficient  to  know  that  such  a  bridge — between  a  precisely  described  %  and  its 

associated  induced  SPT’s — exists.  We  remind  the  reader  that  we  will  address,  in  Section 

V,  this  uncertainty  as  to  which  %  is  actually  executing  during  any  particular  execution  run 

of  application  a  e  ?4. 

11.2  pC  SPT’s 

We  are  now  prepared  to  describe  the  normal  jiiC  execution  of  task  %  as  well  as  the 
disruption  of  the  execution  of  %  by  an  EM  pulse.  In  light  of  our  aforementioned  focus,  we 

will  do  so  from  a  signal-centric  point  of  view.  Consider  then  an  application/system  task  % 

active  during  [0,  TE).  In  what  follows  we  will  also  be  interested  in  some  non-void  “local” 
time  sub-window  [T0,  Tf)  c  [0,  TE)  that  will  designate  the  total  time  extent  of  nEU.  From 

our  signal-centric  point  of  view,  the  task  %  may  be  represented  via  a  |uC  state  that  is 

completely  specified  by  giving,  for  each  relevant  signal  line  of  the  \.iC,  the  SPT  time  history 
on  the  line  during  [0,  TE)  (say  referenced  to  the  line  midpoint  for  definiteness)  hence  also 
during  [T0,  Tf).  That  is,  if  there  are  L  such  signal  lines,  labeled  by  C=  1 ,  ...,  L,  then  the  piC 

state  for  the  time  window  [0,  TE)  is  given  by  a  finite  sequence  <^>[=1  where,  for  normal 
|uC  functionality,  nt  =  ^(rQ:  [0,  TE)  [14,  iow,  W,  high]  is  a  continuous  function  giving  the 
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SPT  time-history  for  signal  line  t,  with  Vn,\ow  representing  “line  low”,  \4,  high  representing 
“line  high”,  and  the  interval  ( \4,  iow,  Vz,high)  accommodating  periods  of  line  transition  from 
low  to  high  or  vice-versa  (and  the  range  of  ?%(%)  has  been  idealized).  For  non-normal  \iC 
functionality,  the  range  of  continuous  function  nt  = 

7re{%-,  7tEU)  on  domain  [0,  TE)  may  extend  beyond  [14,  iow,  V4,  high]  but  is  still  bounded.  In 

general,  it  is  sufficient  to  consider  these  SPT’s — both  normal  and  non-normal — as 
members  of  the  Cartesian  product  XLB{[0,  TE))  =  X[=1  Ce([0,  TE)),  where  Ce([0,  TE))  is  the 
set  of  (real-valued)  continuous  functions  on  [0,  TE)  that  are  uniformly  bounded  there  by  B 
>  0  (i.e.,  [Vfe  Ce([0,  Te))][ ||  f B\),  with  the  specific  value  of  B  being  unimportant  here 
except  that  B  >  max  {|  14,  iow\ ,  1 14,  high |}  (using  “>”  rather  than  “=”  for  the  condition  on  B  in 
order  to  accommodate  non-normal  SPT’s).  Actually,  since  there  are  (almost  certainly) 
<nt>\ =i  e  Xg([0,  Te))  that,  for  physical  reasons,  never  occur  either  as  normal  or  as  non¬ 
normal  jiiC  states,  we  then  denote  by  nL([0,  TE))  ^ XLB{[0 ,  TE ))  the  set  of  all  <^>[=1  e 
Xg([0,  Te))  that  may  in  fact  occur  as  piC  states  (either  normal  or  non-normal).  We  point 
out  that  a  state  </z^(GQ>[=1  for  normal  )liC  functionality  corresponding  to  given  task  %  is 

in  general  not  unique:  Actually,  because  small  but  nonzero  timing  and  amplitude 
variations  in  SPT’s  are  allowed  by  design  without  impacting  normal  jiiC  functionality,  each 

l  =  1,...,  L,  should  be  replaced  with  an  “equivalence  class”  of  such  7re{%)' s  and, 
likewise,  <^(^)>[=1  should  be  replaced  with  an  equivalence  class,  say  &Lnormai[%],  of 
such  <?%(%.)>(  =-\'s  where,  when  <7re{%)>^ =1  then  ~  <^{%)>L i 

iff  they  both  correspond  to  the  successful  jiiC  execution  of  task  %.  (Also,  since  the 

members  of  any  given  normal  class  differ  only  by  intended  “micro”  design  tolerances  and 
not  by  any  “macro”  differences  [e.g.,  totally  differing  line  signal  levels  in  any  given  clock 
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Vfc-period],  we  could  assume  that  <?tLnormal[%]  n  <%Lnormal[%']  =  0  when  support 

the  use  of  the  term  equivalence  class;  for  then  ~  is  a  bona  fide  equivalence  relation  on 
the  set  U*  G  ^  ^Ln0ma/[^]-)  We  wil1’  however,  avoid  equivalence  classes  for  the  most 

part  in  our  main  development  in  Section  III;  rather,  with  any  given  task  %  executing 
normally  we  will  associate  a  unique,  idealized  (as  per  Section  III)  state  representative  of 
its  equivalence  class  <?CLnorma[%],  namely, 

<*?{%)>L  i- 

Consider  next  the  non-normal  states  that  may  occur  when  the  pC  is  subject  to  an  external 
EM  pulse  ;rEM.  To  that  end,  for  pC  task  %  represented  in  normal  execution  by  idealized 

state  <x?{rQ>e=v  let  <?%(%’ xBM)>e=-\  e  IJL{[ 0,  TE))  represent  any  possible  piC  state 
achievable  when  the  piC  is  subjected  to  external  EM  pulse  x™.  In  general,  the 

ttem)>'s  differ  from  <xf{%)>  in  that  the  former  engender  timing  and  amplitude 
variations  of  the  latter.  In  contrast  to  the  normal  case,  there  are  no  compelling  natural 
equivalence  classes  for  the  non-normal  states  <xe{%-,xEM)>^=v  Now  if  state 

<^(rZ^;^-EM)>[=1  £  &iL  normai[% ]  then  this  state  is  said  to  be  deviant  with  respect  to  (wrt) 

task  %  and  written  for  emphasis  as  <^f(<2^;^-EM)>[=1;  further,  this  deviant  state  is  termed, 

in  addition,  disruptive  wrt  task  %  if  it  corresponds  to  mal-execution  of  that  task;  and  it  is 

termed  non-disruptive  wrt  task  %  otherwise.  In  this  connection,  there  might  be  deviant 

states  that — despite  not  being  elements  of  &Lnormai[%] — nevertheless  do  not  correspond 

to  mal-execution  of  task  %  and  are  thus  not  disruptive.  To  capture  the  above 

11 


Approved  for  public  release:  distribution  unlimited. 


considerations  in  our  formulation,  we  define  the  following  subsets  of 

nL^JE))\ecLnormam- 

nvL{%\  [o,rE);  ttem) 

=  {<7re(%-7TEM)>^  e  nL([0,  TE))\&Lnormal[%\  |  <^(<t;;^em)>[=1  is  a  deviant  state  wrt  task  %}, 
nvL-ND(%-[0,TEy,7rEM) 

=  {<**{%■, e  I7*(%-,  [0 ,TE)  ;  nEU)  |  <n*(l Ta',ftEM)>e  =  -\  is  a  non-disruptively  deviant 
state  wrt  task  7}, 

n*c  D(%-’  [0,TeY,xem) 

=  {<nf{%',xEM)>e=i  e  77^^;  [0,T£);  nEU)  |  <^f(rZ^;^EM)>[=1  is  disruptively  deviant  state  wrt 
task  %}, 

for  which  77^^;  [0,  T£);  xEU)  is  the  disjoint  union 

nvL{%-,  [0 ,r£);  XEM)  =  17*' [0 ,r£);  .EV^U;  [0,r£);  7TEM). 

Also,  we  denote  by  nL(%]  [0,  TE))  the  disjoint  union 

nL(%-,  [0,  r£»  =  eiLnor[j.%]  u  "<>(%■,  [o,  r£):  *-EM)  u  n”.  °(<r.;  [o,  te)^u). 

11.3  Probability  of  pC  Upset,  P“pM 

We  may  now  more  carefully  state  in  terms  of  pC  states  (i.e.,  collections  of  SPT’s  on  the 
pC  lines)  what  we  mean  by  pC  upset  probability  for  a  given  %.  Consider  then  the  following 

thought  experiment.  A  fixed  pC  is  executing  a  known,  fixed  % — hence  known  <xf {%)>(= 
1 — and  is  irradiated  by  a  fixed  nEM  acting  over  fixed  time  interval  [7^,  Tf)  c 
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[0,  Te),  with  T0  being  some  fixed  instant  after  %  initiation  (and  with  the  EM  pulse  being 

strong  enough  to  potentially  cause  |uC  upset  but  not  so  strong  as  to  cause  jliC  damage); 
in  this  case  we  also  write  xEU[T0,  Tf\.  Repeated  trials  of  this  procedure  are  conducted  and 

the  outcome  state  <^(%:^BM[T0,  7j])>[=1  e  nL(%\  [0,  TE))  is  observed  and  classified  as 

to  which  of  the  three  disjoint  subsets  of  nL{%\  [0,  TE))  it  belongs.  According  to  our 

comments  in  the  Introduction,  we  take  this  outcome  to  be  probabilistic  (our  EM  pulse 
being  presumed  to  be  in  a  frequency  regime  in  which  this  is  so)  so  that  the  resulting  state 

<^(rZj(;;rEM[7j),  7)])>[=1  and  the  subset  to  which  it  belongs  varies  from  trial  to  trial.  This 
probabilistic  experiment  then  has  outcome  set 

«  -  {ee^VQ,  nvc  ND(%-  [0,  te)\  7tem[T0,  Tf]),  n®L’  [o,  te)-  n™[T0,  t,])} 

and  we  seek  to  theoretically  determine  the  value  of  Prob(77^-  D(cTa,  [0,  TE)\  xEM[T0,  Tf])), 
i.e.,  the  fraction  of  all  trials  in  which  <7re{%-,7rEM[T0,  7j])>[=1  is  disruptively  deviant;  this 
number  will  then  be  declared  to  be  equal  to  pupset  for  task  %  subject  to  nE  M[T0,  Tf]: 

Pupset(%-,xEM[T0,  Tf])  =  Prob(/7®.  D(%:  [0,  Te);  xem[T0,  Tf])). 

We  may  also  think  of  Pupset(%-,xEM[T0,  Tf])  as  (being  the  same  as)  the  probability  of  the 
disruptive  deviation  of  idealized,  normal  SPT  <^7{rQ>(=-[)  subject  to  n  EM[T0,  Tf]: 

Pupset{%:7rEM[T0,  Tf])  =  ProbV'D(<A%)>U  ^EM[4  Tf]). 


In  order  then  to  theoretically  determine  Prob(77^  D(%\  [0,  TE)\  nEU[T0,  Tf])),  we  construct 

a  model  in  which  we  provide,  for  each  possible  idealized,  normal  SPT  <x?{Ta)>c=v  the 

EM-pulse-dependent  probability  of  a  disruptive  deviation  of  <^f(Tjj)>[=1  by  any  given  EM 

pulse  nEM.  Note  that  our  probability  assignment  will  be  a  function  of  only  </rf(Tjj)>[=1  and 
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not  depend  upon  the  resulting  state  <xe(%-,xEM)>(=r  This  approach  frees  us  from  the 

extremely  challenging  (and  indeed  probably  currently  unachievable)  burden  of  precisely 
determining  and  specifying  any  particular  EM-induced  deviant  state  (disruptive  or  not): 

we  need  only  possess  <xf(Ta)>^=)  (which,  if  required  in  detail,  is  in  principle  specifiable — 
but  in  practice  at  best  arduously  so — for  %  from  knowledge  of  the  |uC  design). 

We  will  proceed  by  initially  reducing  consideration  of  the  full  collection  </rf(GQ>[=1  of 
SPT’s  on  all  L  (“important”)  jliC  signal  lines  to  consideration  of  a  single  such  signal  line. 
We  say  that  <^7(rC)>[=1  suffers  a  disruptive  deviation  (i.e.,  that  execution  of  %  is  upset) 

iff  there  is  at  least  one  i  =  1,  ...,  L  such  that  xf{%)  suffers  a  disruptive  deviation  (i.e., 
corresponds  to  mal-execution  of  %)  and  we  then  address  the  determination  of  the 

probability  Prob D(xf{%)-,  xem)  that  any  single  SPT  xf{%)  suffers  a  disruptive  deviation 

induced  by  ^EM.  Of  course,  one  can  then  infer  the  probability  disruptive  deviation  of 

<xf{%)>Le^  by  x™,  i.e.,  Prob^ D(<xf(%)>^=v  xEM),  from 

Prob^  D(<^(rZ^)>[=1;  xEM)  =  1  -  {Prob([;rf  {%)  does  not  suffer  disruptive  deviation 

via  xEU]  &  \xl {%)  does  not  suffer  disruptive  deviation  via  xEU]  ...  &  ...  [xfiT^  does 
not  suffer  disruptive  deviation  via  x™])} 

(which  in  general  is  equal  to  1  -  n[=1  [1  -  Prob^,  D(<^(^)>[=1;  ^EM)]  only  when  the 

individual  single  line  SPT  upsets  are  independent);  but,  because  of  the  need  to  make 
mathematically  precise  the  notional  arguments  just  presented,  as  well  as  to  discuss 

subtleties  that  arise,  we  defer  the  computation  of  such  Pxob^  D{<xf{%)>^=y,  xEU)  to  a 
subsequent  paper.  Thus  the  model  in  this  paper  provides  the  probability  of  disruptive 
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deviation  of  any  particular  idealized,  normal  SPT,  xf(%),  as  induced  by  any  EM  pulse, 
nEM,  acting  over  time  interval  [T0,  Tf)  c  [0,  TE). 

As  a  final  point  in  this  section,  we  remind  the  reader  that  we  have  earlier  briefly  discussed 
the  uncertainty  associated  with  knowing  which  %  is  actually  executing  when  it  is  stated 

that  “application  n  is  executing”  and  have  alluded  to  treating  that  stochastically  (in  Section 
V).  That  stochasticity  is,  however,  completely  different  from  that  discussed  in  the  previous 
paragraph  in  which  the  actually  executing  %  is  taken  to  be  fully  known  but  in  which  the 

upset  of  7r?{%)  is  stochastic. 

III.  SPT  (a71)  and  EM  Pulse  (;rEM)  Specification  and  Interaction 

As  indicated  in  the  previous  section,  we  are  interested  in  determining 

Prob  ;  xEM)  when  %  is  presumed  known,  hence  (<^f(T^)>[=1  being  known  as 

well.  We  will,  however,  more  generally,  determine  Probp  (^;  xEM)  for  every  possible 

idealized,  normal  SPT  n71  in  [0,  TE).  Since  for  l  =  1 ,  2,  ...,  L,  xf(%)  is  one  of  these  n71' s 

then  Prob^ D(;rf(%)m,  xEU)  will  follow  from  our  more  general  results.  To  this  end,  we 

carefully  specify  in  this  section  our  idealized  SPT’s,  n 71 ,  and  our  EM  pulses,  xEM,  as  well 

as  characterize  the  overlap  in  time  between  them.  We  consider  a  pC  with  a  fundamental 
clock  period  of  xc  >  0  and  denote  xc  1/2  =  xc/2;  the  latter  specifies  our  basic  time  scale.  We 
thus  consider  the  global  time  interval  [0,  TE)  to  be  partitioned  into  N  >  1  basic  time 
segments  of  length  xc  1/2  according  to 

[0,  T E)  =  U^o  [nx^/2,  {n  +  1)xc1/2),  (1) 

where  we  stipulate  without  practical  consequence  that  TE  is  an  exact  multiple  of  xc  1/2;  i.e., 
Te  =  A/xc  1/2.  We  will  specify  our  SPT,  n71,  in  terms  of  the  partition’s  basic  time  segments. 
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111.1.  Individual  SPT  Elements 

We  consider  any  idealized  SPT,  n71  in  [0,  TE),  to  be  constructible — in  a  graph  of  signal 

magnitude  vs.  time — as  a  concatenation  of  shortly-to-be-specified  basic  elements.  Our 
basic  elements  are  of  two  distinct  classes,  which  two  classes  we  term  as  “primitive  SPT 
elements”  and  as  “Vz-period-segmented  SPT  elements”.  There  are  four  differing 
elements — ’’types” — in  the  former  class  and  eight  differing  element  types  in  the  latter  class 
and  in  this  section  we  precisely  specify  all  these. 

Consider  firstly  the  primitive  SPT  elements.  They  are  of  the  following  types:  (A)  signal- 
high,  having  constant  value  -oo  <  VK,  high  <  oo  during  some  positive  length  time  interval 
(where  14,  high  is  frequently  taken  to  be  positive);  (B)  signal-low,  having  constant  value  -oo 
<  V4,  low  <  V4,  high  during  some  positive  length  time  interval  (where  14,  iow  is  frequently  be 
taken  to  be  zero);  (C)  signal-low-to-high-transition,  a  straight  line  from  termination  of  a 
signal-low  element  to  commencement — at  some  positive  time  beyond  that  termination — 
of  the  very  next-in-time  signal-high  element,  and  with  that  straight  line  having  finite, 
positive  slope  [( 14,  high  -  W,/ow)/2)]/aTc1/2,  with  a  >  0;  and  (D)  signal-high-to-low-transition, 
a  straight  line  from  termination  of  a  signal-high  element  to  commencement — at  some 
positive  time  beyond  that  termination — of  the  very  next-in-time  signal-low  element,  and 
with  that  straight  line  having  finite,  negative  slope 

-[(VSr.  high  Mr,  /ow)/2)]/|3tc  1/2,  with  p  >  0.  In  practice,  it  usually  occurs  that  a,  p  «  1 ;  less 
stringently  we  stipulate  only  that  0  <  a  +  p  <  1  (so  that  0  <  a  <  1  and  0  <  p  <  1 ).  Examples 
of  primitive  SPT  elements  are  illustrated  in  the  SPT  of  Fig.  1.  Note  that  primitive  SPT 
elements  of  Types  A  and  B  may  (partially)  occupy  more  than  a  single  basic  time  segment 
while  elements  of  Type  C  and  D  always  (partially)  occupy  more  than  a  single  basic  time 
segment. 

Consider  next  the  V^-period-segmented  SPT  elements.  Now  in  general,  the  portion  of  any 
SPT  occurring  during  any  single  basic  time  segment  [mc1/2,  (. n  +  1)tcJ/2)  is  constructible 
by  concatenation  of  portions  of  several  of  the  above  primitive  SPT  element  types.  We 
illustrate  in  Fig.  2  all  eight  possibilities  (types)  for  the  structure  of  an  SPT  during  (without 
loss  of  generality)  the  zeroth  basic  time  segment  {n  =  0  in  Eq.  (1 ));  and  we  refer  to  these 
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eight  as  “Vfe -period-segmented  SPT  elements".  As  an  example  of  the  constructability  of 
a  1/2-period-segmented  SPT  element  from  primitive  SPT  elements,  note  that  the  !/2-period- 
segmented  SPT  element  of  Type-4  may  be  regarded  as  the  concatenation  of  a  primitive 
element  of  Type-A  with  half  of  each  of  primitive  elements  of  Type-C  and  Type-D. 

The  reason  we  have  introduced  both  primitive  elements  and  V^-period-segmented 
elements  is  as  follows.  While  the  latter  are  more  natural  for  describing  an  SPT  in  time,  it 
nevertheless  seems  to  us  that  the  former  are  more  natural  for  assigning  disruptive 
deviation  probabilities.  Our  ultimate  method  of  assigning  disruptive  deviation  probabilities 
to  single  V^-period-segmented  SPT  elements  (and  therefrom  to  an  /V-element  SPT)  will 
be  to  first  assign  disruptive  deviation  probabilities  to  basic  (see  Section  IV.2.1)  primitive 
SPT  elements  and  therefrom  compute  the  disruptive  deviation  probabilities  for  the  V2- 
period-segmented  SPT  elements  (Section  IV. 2).  However,  in  what  follows  we  will  also 
initially  take  the  simpler  approach  of  assigning  disruptive  deviation  probabilities  directly 
to  the  !/2-period-segmented  SPT  elements  themselves  (Section  IV.  1). 

To  proceed,  we  label  the  eight  Vb-period-segmented  SPT  elements  in  Fig.  2  by  sy, 
j=  1,...,  8,  as  indicated  specifically  in  that  figure.  In  fact,  each  such  sy  is  a  mapping  from 
[0,  tCii /2)  into  [Vn,  low,  Vn,  high]-  Explicitly,  denoting  A 14  =  Vn,  high  -  14,  ,ow  >  0  and 
Vn,  mid=  {Vn,  high  +  Vn,  i0W)I2  we  have  the  following: 

Sl(0  =  Vn,  high  t  e  [0,  Xc1/2)  (2) 


s2{t)  =< 


(Al4/axc)  t  +  Vn,  mid  if  t  e  [0,  axCj1/2) 

Vn,  high  if  t  e  [axc  i/2,  XCii/2)  , 


Vn. 


high 


s3(0  =  i 


—  (Al4/pXc)  t  +  (A14/2P)  +  Vn,  mid 
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if  t  e  [0,  (1  -  p)xCi1/2) 

if  te[{  1  -  P)xc  i/2,  xCj1/2), 


(3) 


(4) 
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r 

(A  14/ CtXc)  t  +  V4,  mid 


s4{t)  =  J 


Vn 


x,  high 


-  (Al4/pxc)  t  +  (A  V4/2(3)  +  \4, 

mid 


if  f  e  [0,  axCj1/2) 

if  f  e  [axc j/2,  (1  -  P)xc  1/2)  (5) 

if  f  e  [(1  -  P)xc  i/2,  xc 

,1/2) . 


55(f)  =  14,  fow  f e  [0,  "4,1/2) 


(6) 


se(0  = 


—  (Al4/pxc)  t  +  Vjt,  mid  if  f  e  [0,  Pxc  i/2) 

V4  ,  low  if  t  e  [pxCii/2,  xc  i/2)  , 


(7) 


s7(0  = 


V4,  fow  if  f  e  [0,  (1  -  a)xc  1/2) 

.  (Al4/axc)  t  —  (Al4/2a)  +  14,  mid  if  t  e  [(1  —  a)xc  i/2,  xc  i/2). 


(8) 


ss(0 


-  (A\4/pxc)f+  \4 


rr,  mid 


=  < 


V, 


x,  low 


V 


(A  \4/axc)  f  -  (A  14 /2a)  +  \4 


x,  mid 


if  t  e  [0,  pxc  i/2) 


if  t  e  [pxc>1/2,  (1  -  a)xc  1/2) 
if  f  e  [(1  -  a)xc  i/2,  xc  1/2). 


0) 


It  is  convenient  to  extend  the  domain  of  each  of  these  V^-period-segmented  SPT  elements 
to  all  of  [0,  Te)\  we  take 


Sj{t)  if  t  e  [0,  xc  1/2) 


s?(t)  -  A 


(j=  1 . 8) 


(10) 


0  if  t  e  [0,  Te)  \  [0,  xc  1/2) 

and  refer  to  these  as  well  (with  no  confusion  generally  resulting  from  doing  so)  as 
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1/2-period-segmented  SPT  elements.  (The  possibility  that  VK,i0W  or  VK,  high  may  have  value 
zero  is  purely  coincidental  in  relation  to  the  possible  value  0  taken  on  by  sy°(f)  in 
[0,  Te)  \  [0,  tc  1/2).)  Further,  to  describe  the  SPT  in  basic  time  segments  other  than  the 
zeroth  one,  we  define,  for  n  =  1 N-  1  and  t  e  [0,  TE),  the  nth-time-segment  V2-period- 
segmented  SPT  elements 


5/(0=  i 


s°(f-  nxc1/2)  if  t  e  [nxc1/2,  (n  +  1)  xc1/2) 

0  if  t  e  [0,  Te)  \  [nxc1/2,  (n  +  1  )xc  1/2); 


(7=1.-.  8)  (11) 


and,  for  each  n  =  0,...,  A/ -  1,  we  denote  the  set  of  all  such  V^-period-segmented  SPT 
elements  as 


(12) 


111.2.  Consecutive  V2-Period-Seqmented  SPT  Elements:  Arbitrary  Idealized  SPT’s 

We  now  use  the  individual  V^-period-segmented  SPT  elements  of  the  previous  subsection 
to  construct  an  arbitrary  A/-element  SPT  having  domain  [0,  TE)  partitioned  into  basic  time 
segments  as  per  Eq.  (1),  with  this  /V-e!ement  SPT  consisting  of  a  sequence  of  N 
consecutive-in-time  single  SPT  elements,  where  the  nth  such  single  element  [n  =  0,...,  N 

-  1)  of  the  sequence  belongs  to  S"  of  Eq.  (12).  In  Fig.  3,  we  show  an  8-element  SPT 
(which,  incidentally  consists  of  all  eight  V^-period-segmented  SPT  element  types).  In 
general,  any  such  consecutive-A/-element  SPT  can  be  written  as  a  sum,  over  n  =  0,...,  N 

-  1 ,  of  individual  SPT  elements  from  the  S^'s;  to  wit,  as  a  sum  Z^fo  sfn  where  s f  e  Sj1 

with  jn  e  {1 ,  ...,  8}.  On  the  other  hand,  it  is  not  the  case  that  every  sum  of  this  form  is  an 
SPT:  While  in  any  given  basic  time  segment  (indexed  by  n)  each  of  the  eight  V2-period- 
segmented  elements  (sj1  for yn  =  1 ,...,  8)  is  possible,  nevertheless  the  possible  pairs,  <sfn , 

s?n\\>,  of  V^-period-segmented  elements  occurring  in  consecutive  time  segments  are 
constrained  by  the  requirement  that  the  SPT  be  a  continuous  function  in  [mc  1/2,  (n  + 
2)xc,i/2)  and  differentiable  at  t  =  {n  +  1)xc1/2  (but  clearly  need  not  be  differentiable  in  all  of 
[mc  i/2,  {n  +  2)xc1/2)).  The  required  differentiability  at  t  =  {n  +  1)xc1/2  serves  to  exclude  any 

19 


Approved  for  public  release:  distribution  unlimited. 


pair  whose  sum  is  continuous  yet  does  not  occur  as  part  of  any  SPT,  e.g.,  type-2  following 
type-3,  or  type-2  following  type-4,  or  type-4  following  type-3.  It  is  straightforward  to  verify 
(from  Fig.  2,  for  example)  that  given  any  n  =  0,...,  N-  1  and  any  s f  e  S"  {jn  e  {1 ,  ...,  8}) 

there  are  then  precisely  two  members  of  S1n+ 1  for  which  the  “two-element”  SPT 
[ s fn  +  s?n\\]  \[nxCV2,  {n  +  2)xc  1/2)  is  admissible  in  the  above  sense.  For  example,  if  jn  =  1 
then  we  must  have  jn+1  e{1 , 3}.  Furthermore,  this  relationship  is  independent  of  n:  if 

n'  *  n  and  }n,  =  1  then  jn,  + 1  e  {1 , 3};  and  this  n-independence  holds  as  well  for 

jn  =  2,...,8.  Thus,  in  general,  for  each  j=  1 ,  ...,  8  we  may  write  (j)+  and  (j) \  for  those  two 

indices  for  which — and  only  for  which — the  “two-element”  SPT’s 

[Sj  +  s^]  ([0,  2xc1/2)  and  [s°  +  s^+]  ([0,  2xCi1/2)  are  admissible  (hence  so  are 

[Sjn  +  S(jy,")+]  f[nxCj1/2,  [n  +  2)xc  1/2)  and  [sj7,  +  s^l]  f[nxCj1/2,  [n  +  2)xc>1/2)  admissible  for  all 
n  =  0,  ...,  N-  1).  We  also  write 

^  ={<} U*>  |  7=1,...,  8  and  J+ =  {j)t  ory+=(y)J}.  (13) 

We  may  now  characterize  all  idealized,  “normally  executing”  /V-element  SPT’s  in  [0,  TE) 
as  the  members  of  the  set 

J7110,  X£]  =  {**|**  =  2;  s?n  with  ( Vn  =  0, . . N  -  1  )(<  y„,  j„  t ,  >  <=  J* ))  (14) 

where 

y  [mc>1/2,  (n  +  1)xc1/2)  — »  [14,  /ow,  ^4,  /i/g/?] ■  (15) 

The  members  of  77^[0,  TE]  may  each  have  one  of  only  eight  possible  individual  V2-period- 

segmented  SPT  elements  occupying  [0,  xc  1/2)  and  for  each  of  these  eight  there  are  two 
possible  individual  V^-period-segmented  SPT  elements  occupying  each  of 

[nxc  1/2,  {n  +  1)xc1/2),  n=  1,...,  N-  1;  hence  #(77^)  =  8*2/v~1  =2W+2.  Further,  since  each 
e77^[0,  Te]  is  completely  characterized  by  its  associated  sequence 
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<j0,  jv...,  jN-r>  (when  a,  (3,  V4,  iow,  and  14,  high,  are  fixed)  then,  when  necessary,  we  will 
denote  that  SPT  by  n*[j 0,  /, . 


III. 3.  EM  Pulses 

An  EM  pulse,  ^EM,  will  be  taken  to  be  a  rectangular-envelope-modulated  monochromatic 
sinusoidal  wave  of  frequency  co  >  0,  having  amplitude  \/EM  >  0  during  time  interval  [7^,  Tf) 

c  [0,  Te)  (with  Tf>T0)  and  amplitude  zero  during  [0,  TE)  \  [T0,  Tf).  Explicitly, 


7TEM{t)  = 


Vem  sin[o(f-  70)] 
0 


iffe[70,  Tf) 


if  f e  [0,  Te)\[T0,  Tf), 


(16) 


where  we  stipulate  without  practical  consequence  that  Tf  -  T0  is  an  integral  multiple  of  tc/co 
so  that  nEU  is  continuous  at  both  T0  and  Tf  (we  attribute  no  significance  to  whether  such 
multiple  is  odd  or  even).  We  will  also  make  use  of  the  upper-half  envelope  of  nEU, 
namely, 


t.em(0  =  i 


^em  t  e  [T0,  Tf) 

0  if  f  e  [0,  Te)\[T0,  Tf). 


(17) 


Further,  when  necessary,  we  will  write  both  and  ^+EM  more  explicitly  as,  respectively, 
kem[Tq,  Tf,  \/EM,  co]  and  ;r+EM[70,  T„  \/EM]. 


111.4.  Overlap  in  Time  of  ^and  nEU 

Given  x™[Tq,  Tf,  \/EM],  we  will  be  interested,  for  use  in  the  remainder  of  the  paper,  in 
characterizing  the  temporal  overlap  of  ;r+EM[7^,  Tf,  14m],  when  Tf  >  T0,  and  any  n71  e 

ZT^[0,  Te],  which  overlap  we  indicate  symbolically  also  as  n71  n  xEU\  i.e.,  n71  n  xEU  = 

[T0,  Tf).  To  this  end,  we  define 
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and 


nmin[T0,  Tf]  =  min{  n  e  {0,  N-  1}  |  [nxc>1/2>  {n  +  1)  xc,1/2)  n 


^maJ7o-  s  max{  n  e  {0,  ...,  A/  -  1}  |  [nxCi1/2,  (n  +  1)  xc>1/2)  n 


where,  of  course,  N-  1  >  n  >  n->  0.  Then 

IllaX  min 


[To,  Tf)  -  { [nmjnzcV2,  (nmjn+  1)xc,i/2)  n  [To,  Tf ) } 


u  Ur^'i  [mc,1/2,(n+1)Tc,1/2) 


[T0,  Tf)  ^0} 
[T0,  Tf)  ^0) 


(18) 

(19) 


(20) 


U  (  [nmaxXcM2'  (nmax+  ”')xc,1/2)  n  [T),  Tf)  }, 

since  [nxci/2,  (n  +  1  )xc  1/2)  [T0,  Tf)  =  [mcV2,  ( n  +  1)xc1/2)  for  n  =  nmjn  +  1,...,  nmax~  1- 

However,  in  general,  [nminx Cj1/2,  (nm//J+  1)xc,1/2)  n  [70,  7f)  ^  [70,  7f)  and 

[nmaxTc,i/2.  (nmax+  1)xc,i/2)  n  [70,  7f)  ^  [70,  Tf).  In  other  words,  the  overlap  in  time  of 
x™[T0,  Tf,  \/EM]  and  any  7r^  s  77^(0,  Te]  is  specified  by  the  collection  of  nmax  —  ^rnin  ~  ^ 

“complete”  time  segments  contained  inside  the  middle  (“big  union”)  term  in  Eq.  (20) 
together  with  the  two  generally  “partial”  time  sub-segments  given  by  the  first  and  last 
terms  in  that  equation.  Our  interest  in  so  specifying  this  overlap  derives  from  the 
circumstance  that  in  the  next  section  we  will  ascribe  EM-pulse  disruption  probabilities  to 
each  of  the  eight  possible  V^-period-segmented  SPT  elements  and  therefrom  infer  an  EM- 
pulse  disruption  probability  for  any  entire  /V-element  SPT  (i.e.,  any  x*  e  77*  [0,  Te]) 

constructed  from  those  V^-period-segmented  SPT  elements.  Each  such  V2-period- 
segmented  SPT  element  of  course  temporally  “occupies”  a  complete  basic  time  segment; 
but,  because  of  the  possibility  that  the  leading  and  trailing  ends  of  7T71  n  ;xEM  each  only 

partially  occupy  a  basic  time  segment,  we  must  in  addition  ascribe  EM-pulse  disruption 
probabilities  to  all  possible  partial  V^-period-segmented  SPT  elements  formable  from  any 
of  the  eight  complete  V^-period-segmented  SPT  elements.  In  this  connection,  it  will  also 
be  useful  to  denote 


Ik71  n  ;rEM[  =  [nmjnzcV2,  {nmax+  1  )xc,i/2)> 


(21) 
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this  quantity  being  the  minimal  extension  of  the  interval  [T0,  Tf)  to  a  union  of  an  integer 
number  of  complete  basic  time  segments,  where  the  interval  [T0,  Tf)  =  n71  n  x™  is  not  in 
general  such  a  union  but  rather  properly  contained  inside  ]^n  ^EM[ . 

IV.  Probability  of  a  ;rEM-lnduced  Disruptive  Deviation  of 

In  this  section  we  formulate  our  expression  for  the  probability  of  a  ^EM-induced  disruptive 
deviation  of  any  single  signal  line  carrying  any  of  the  possible  SPT’s,  n71  e 

77^[0,  Te\.  We  firstly  consider  disruptive  deviation  probability  assignments  to  each  of  the 

individual  V^-period-segmented  SPT  elements  sj7  and  do  so  from  two  points  of  view:  (1) 
direct  assignment  to  each  of  the  eight  V^-period-segmented  SPT  elements  and  (2)  indirect 
assignment  to  each  of  the  eight  V^-period-segmented  SPT  elements  via  direct  assignment 
of  disruptive  deviation  probabilities  to  each  of  the  four  primitive  SPT  elements.  We  then 
consider,  secondly,  disruption  probability  assignment  to  a  single,  arbitrary  /V-element  SPT 
built  from  the  sj1' s.  The  connection  between  these  two  entities  is  this:  Given  an  SPT 
composed  of  a  specified  set  of  V^-period-segmented  SPT  elements  then  we  say  that  the 
SPT  suffers  a  disruptive  deviation  iff  at  least  one  of  its  elements  suffers  a  disruptive 
deviation.  Thirdly,  we  consider  (in  Section  V)  the  disruptive  deviation  probability 
assignment  to  a  single  signal  line  in  case  the  actual  SPT  on  the  line — being  an  SPT  from 
the  collection,  njj,  of  all  possible  A/-element  SPT’s  that  may  be  carried  on  the  line — is 

known  only  stochastically. 


IV.1 .  Individual  V^-Period-Seqmented  SPT  Elements:  Direct  Probability  Assignment 
IV.1 .1 .  Full  Overlap  of  nEM  and  s" 

Consider  firstly  any  EM  pulse  ;xEM[0,  Tf,  VEM,  co]  such  that  [T0,  Tf )  n  [0,xc1/2)  =  [0,xc1/2),  so 
that  [T0,  7j)  3  [0,xc  i/2)  (and  Tf  >  xc1/2),  and  further  consider  l/2-period-segmented  SPT 
elements  sf  e  S°,  j  =  1,...,  8.  In  this  case  we  denote  by  p°  fu//[V^M,  co]  €  [0,  1]  the 
probability  that  element  sy°  suffers  a  disruptive  deviation  induced  by  ^EM.  (We  have 
suppressed  the  implicit  dependence  of  p°  M  upon  xc1/2,  V4,  iow,  V4,  high,  a,  and  p.)  Note 
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that  we  are  considering  this  probability  to  be  independent  of  Tf  as  long  as  Tf  >  xCj1/2:  The 
differences  in  the  probabilities  of  ^EM-induced  disruptive  deviations  of  multi-segment 
SPT's  for  differing  such  ^EM 's — the  ^EM 's  differences  arising  from  the  differences  in  their 
Tf' s  (when  all  of  their  T/s  are  strictly  greater  than  xc1/2) — will  be  manifested  in  our  model 
in  the  V^-period-segmented  SPT  elements  beyond  n  =  0. 

Next,  for  each  n  =  1 N-  1 ,  consider  any  EM  pulse  ^EM[7^,  Tf,  \/EM,  co]  such  that 
[T0,  Tf)  n  [nxc1/2,  (n  +  1)xc>1/2)  =  [nxc1/2,  ( n  +  1)xci/2),  so  that  [7^,  7))  z? 

[nx Cji/2,  ( n  +  1)xc1/2)  (and  T0  <  nxc  1/2  and  Tf>{n  +  1)xc1/2),  and  further  consider  V2-period- 
segmented  SPT  elements  sj1  e  S°,j=  1,...,  8.  In  this  case,  we  denote  by  p" fut{  VEM,  co]  the 
probability  that  element  sj7  suffers  a  disruptive  deviation  induced  by  ^EM  where,  once 
again,  we  are  considering  this  probability  to  be  independent  ofT0  and  Tf  as  long  as  T0< 
nxc  i/2  and  Tf>{n  +  1)xc1/2:  The  differences  in  the  probabilities  of  xEU  disruptions  of  multi¬ 
segment  SPT's  for  differing  such  ^EM 's  will  be  manifested  in  our  model  in  the  V2-period- 
segmented  SPT  elements  previous  to  and  subsequent  to  the  nth.  Further,  we  assume  that 
Pi,  fuiiiVEw  «]  =  Pi,  fuiil ^EM’  “L  7=1 . 8;n=1 . A/-1;  (22) 

i.e.,  for  each  one  of  the  eight  V^period-segmented  elements,  its  probability  of  suffering  a 
disruptive  deviation  is  independent  of  where  that  element  occurs  in  the  N-element  SPT, 
as  long  as  the  element  is  completely  encompassed  by  nEM,  this  probability  depending 
rather  only  upon  which  one  of  the  possible  eight  the  element  actually  is.  We  term  this 
property  “ modularity it  appears  to  us  to  be  a  necessary  assumption  if  we  are  to  model 
the  upset  of  all  members  of  the  collection  of  all  possible  SPT’s  in  [0,  TE)  without  treating 
each  such  member  as  an  individual  having  no  commonality  with  any  other  member  of  the 
collection  with  respect  to  response  to  an  EM  excitation. 

IV.  1 .2.  Partial  Overlap  of  ;rEM  and  sj1 

Suppose  next,  contrary  to  the  first  paragraph  of  Section  IV.1 .1 ,  that  ^EM[7^,  Tf,  VEM,  co]  is 
such  that  [T0,  Tf)  n  [0,xc1/2)  ^  [0,xcJ/2)  and  [T0,  Tf)  n  [0,  xcJ/2)  t  0,  with 

24 


Approved  for  public  release:  distribution  unlimited. 


VT0’  Tf)  =  Length {[T0,  Tf)  n  [0,  tCj1/2)}/xCi1/2; 


(23) 


note  0  <Y0(T0,  Tf)  <  1 .  In  this  case,  we  assume  that  the  probability  that  element  sy°  suffers 
a  disruptive  deviation  induced  by  nEU  is  given  by 

Pf.  partial  t  ^EM’  03  ’  ^0>  =  {/(V^O’  Tf),  pfjun[VEM,  Co])  pfjun[VEM,  Co],  (24) 

where  g\  [0,  I]2  ->  [0,  1]  is  some  yet-to-be-specified  function  of  Y0  and  Py°,  full  that  is  non¬ 
decreasing  in  Y0  for  each  possible  fixed  p°  fc//  e  [0,  1]  and  which  further  satisfies 

(Vp“  „  e  [0, 1])<g(0,  Pi  m)  =  0  and  g(1 ,  p“  full)  =  1 ).  (25) 

This  non-decreasing  behavior  of  g  in  Y0  seems  reasonable  to  us  but  may  in  fact  not  obtain 

in  reality;  nevertheless,  we  adopt  this  behavior  in  our  present  model  until  experiment 
dictates  otherwise.  Note  that  Eq.  (24)  stipulates  that — as  far  as  the  interval  dependency 
of  g  is  concerned — it  is  only  the  length  of  the  interval  [7^,  Tf)  n  [0,  tCj1/2)  that  determines 

Py°,  partial ancl  not  the  specific  location  of  this  interval  inside  of  [0,  xc1/2).  (We  will  modify  this 
assumption  later  on,  in  Section  IV. 2 — for  the  present  it  will  suffice  for  our  purpose.)  Also 
note  that  the  only  dependency  of  g  upon  j  is  via  p°  m — the  functional  form  of  g  is  the 
same  for  each  j  =  1 ,  ...,  8  (i.e.,  g  does  not  carry  a  subscript  “j").  Further,  the  inclusion  of 
values  0  and  1  for  Y0,  in  addition  to  those  allowed  values  of  Y0  specified  immediately 
following  Eq.  (23),  will  allow  us  to  incorporate  the  limiting  cases  of  [7^,  Tf)  n  [0,  tc>i/2)  =  0 
and  [7^,  Tf)  n  [0,  tc1/2)  =  [0,  tc1/2),  which  cases  we  will  employ  shortly  in  explicitly 
determining  g. 

Similarly,  for  each  n=  1 ,...,  N -  1 ,  consider  any  EM  pulse  7rEM[T0,  Tf,  \/EM,  co]  such  that  [7^, 

Tf)  n  [nxc1/2,  ( n  +  1)xc1/2)  ^  [r?xc,i/2>  in  +  "I  )xc,i/2)  ar|d  [To,  Tf)  n  [nxCj1/2,  ( n  +  1)xCj1/2)  ^  0, 

with 

yn{T0,  Tf)  =  Length {[T0,  Tf)  n  [nxCj1/2l  (n  +  1  )xc>1/2)  }/xc,1/2;  (26) 
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0<Y n(T0,  Tf)  <  1  (n  =  1 N-  1). 


(27) 


Then  we  assume  that  the  probability  that  element  s"  suffers  a  disruptive  deviation  induced 
by  7rEM  is  given  by 


Pj,  partial  i  ^EM’  Tqj  7)]  ~  ^0n(7 0>  7)),  Pj,  fulli  ^EM’  ®])  Pj,  fulli^EW  ®]) 


-  o,  7)),  Py°  ftj//[^EM>  0)])  Py°  fty//[^EM>  0)] 


(28) 


where  we  have  used  here  one  and  the  same  function  ^  for  all  n  =  1 N-  1 ,  that  function 
also  being  the  one  that  we  used  for  n  =  0.  We  emphasize  that  we  have  thus  taken  the 
functional  form  of  g  to  be  independent  of  both  j  and  n.  Also,  as  in  the  n  =  0  case,  we 

allow  in  addition  ln  =  0,  1  . 

We  now  address  the  determination  of  the  function  g.  For  the  purpose  of  this  discussion, 
we  will  suppress  the  subscript  “j”  and  superscript  “n”  on  p-  full  and  p"  part/a/ ,  as  well  as  the 
subscript  “n”  on  Jn{T0,  Tf),  and  also  suppress  all  their  arguments;  we  merely  write  pfull , 

Ppmiah  ancl  7  for  these,  with  n  and  y  understood  to  be  fixed  but  arbitrary  in  their  allowed 
ranges.  Also,  for  proper  non-void  subinterval  /  =  [7^,  Tf)  n  [mc1/2,  (n  +  1)tc1/2)  of  basic 
time  segment  [nxc1/2,  (n  +  1)xCj1/2) — with  the  former  being  a  subinterval  whose  p^artiai  is 

that  of  interest  and  whose  length-fraction  is  Y,  and  with  the  latter  having  disruption 
probability  pfull  when  fully  encompassed  by  nEU — we  write 

PpanJl’  Pj  =  PJPm-  (29) 

(Recall,  as  indicated  earlier,  that  two  different  subintervals  inside  the  same  basic  time 
segment  [so  that  both  of  these  subintervals  are  referenced  to  the  same  pfull]  and  having 

the  same  Y  also  have  the  same  p^artiai{y,  pfull).)  Further,  for  r  =  [nxc1/2,  {n  +1  )xc  1/2)  \  I,  with 

having  length-fraction  1  -  Y>  we  take  its  P  jal  as 

<3°) 


26 


Approved  for  public  release:  distribution  unlimited. 


despite  the  fact  that  I~  may  sometimes  not  be  an  interval  but  rather  the  union  of  two 
intervals — we  still  assign  a  Ppartial  to  7~  based  solely  upon  its  total  length.  This  is  in  keeping 
with  our  implied  basic  assumption  that  Ppartja,  for  a  (measurable)  subset  of  a  given  basic 
time  segment  (with  the  latter’s  pfull  hence  also  given)  depends  solely  upon  the  length  of 
that  subset — and  nothing  else. 

To  proceed  with  the  determination  of  g,  consider  a  random  experiment  in  which 
V^-period-segmented  SPT  element  s f — which  element  when  fully  encompassed  by 
nEM[T0,  Tf]  (suppressing  VEM  and  co  here  as  well)  exhibits  upset  probability  pfull — is  in  fact 
irradiated  by  EM  pulse  ttem[T0,  Tf]  having  [T0,  Tf)  =  [nmax tCj1/2,  ( nmax+  1)xc,1/2).  With  some 
specified  I  of  interest  as  per  above  and,  in  particular,  having  length  fraction  y,  one 
observes  and  records  the  subset — I  or  7~ — during  which  a  disruptive  deviation  of  the  SPT 
element  s f  first  takes  place,  making  allowance  as  well  for  the  additional  possibility  that  no 
disruptive  deviation  of  the  SPT  element  sf  takes  place  at  all — which  outcome  is  signified 
by  W  (i.e.,  no  disruptive  deviation).  For  this  experiment,  the  outcome  set  is 

Q  =  {/,  7~,  W}  with  probability  measure  (on  2n — the  power  set  of  Q)  specified  by  P({7})  = 
Pn  P({1  ~})  =  Pr ,  and  P({W})  =  1  -  {Pj  +  pr).  Now  in  terms  of  the  quantities  Ppartiai  and 
pfull  above  we  must  have 


Pz  +  (1  -  P/)  Pr 

II 

II 

! 

+ 

1 

l 

(31a) 

Pl=  Ppartial^’  Pful) 

and  A-  =  Pp,«<1“y'PM>- 

(31b) 

The  first  equality  in  Eq.  (31a)  reflects  the  following:  The  fraction  of  experimental  trials  in 
which  a  disruption  takes  place  in  s" — namely  pfu—  is  equal  to  the  fraction  of  trials  in  which 
the  outcome  is  7—  namely  — plus  the  fraction  of  trials  in  which  the  outcome  is  not  7  but 

is  r — namely  (1  -  p,)pr ;  further,  as  expressed  by  the  second  equality  in  Eq.  (31a),  this 

sum  is  invariant  under  the  interchange  of  I  and  7~.  We  will  use  this  reasoning  again  in 
the  sequel.  Continuing,  we  use  Eqs.  (29)  and  (30)  in  conjunction  with  Eqs.  (31)  to  find 
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<S(t.  Pj  +5(1-1-.  pj  -  Pu,  fffir.  Pj  5(1  -  y.  Pj =  1  ■  (32) 

Abbreviating  p=  pfull,  we  thus  seek  solutions  Cp  [0,  I]2  ->  [0,  1]  of 

£(Y>  P)  +  £(1  ~y>p)~  p£( T»  p)£(1  -  V.  P)  =1  (33) 

for  all  y,  p  e  [0,  1],  subject  to  the  conditions  of  Eq.  (25),  namely, 

(Vp  e  [0,  1  ])(5(0,  p)  =  0  and  g(1 ,  p)  =  1 ).  (34) 

Note  that  Eq.  (33)  demands  that 

pgHV2,p)-2g(V2,p)  +  i  =o  (35) 


for  all  p  e  [0,  1]  and  that  this  equation  has  one  and  only  one  of  its  solutions  for  g  having 
value  in  [0,1],  namely  the  solution 


[i  -(i  -P)m]/p 


Vz 


if  P  e  (0,  1] 
if  p  =  0. 


(36) 


We  now  exhibit  a  solution  to  Eqs.  (33)  +  (34).  (This  solution  is  not  unique — this  will 
become  evident  shortly.)  We  have 


[1  -  (1  -p)5]/p  if  p  e  (0,  1 )  &  8  e  [0,  1] 


£(S,p)=  *1 


8 

0 


if  p  =  0  &  8  e  [0,  1] 
if  p  =  1  &  8  =  0 


1 


\\p=  1  &  8  e  (0,  1] 


(37) 


where  we  substitute  S  =  YorS=1  -  Yas  applicable.  (Eq.  (37)  may  be  written  more 
succinctly  but  for  clarity  we  have  chosen  not  to  do  so.)  Using  Eq.  (34)  along  with  the 
observation  that  (dg/d8)(8,  p)  >  0  for  p  e  (0,  1)  and  8  e  [0,  1],  we  see  that  indeed  £(8,  p) 
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e  [0,  1]  in  [0,  I]2.  Further,  g(8,  p)  is  continuous  in  [0,  1]2\{<0,  1>}  in  each  of  its  variables 
(we  use  the  notation  ( ,  >  for  ordered  pair)  since  for  all  8  e  [0,  1]  we  have 
limJp-H)+  {[1  -  C1  -pf]/p}  =  5  =  S(8>  0)  and 


lim 


p^>  i 


{[i  -  (i  -p)6]/p}  =  A 


0  if  8  =  0 
1  if  8  e  (0,  1] 


=  #8,1) 


(38) 


while,  on  the  other  hand,  lim5^0+  #8,  1)  =  1  *  #0,  1).  The  expression  for /^art(a/  resulting 
via  Eq.  (29)  from  g  given  in  Eq.  (37),  this  former  equation  now  being  properly  and 
completely  rewritten  (in  view  of  Eq.  (22),  hence  Eq.  (28))  as 


P"  partial  n(T0’  Tf )>  P full )  =  Tf )>  Pj,  full )  P\ 


y,  full  » 


(39) 


is 


Pj,  partial^ ni^O,  Tf) ,  Pj  ^//)  —  "S 


(i  -  pI  m)1 n 

0 

1 


if  ftiff  e  [0,  1)  &  yne[0,  1] 

if  Py°,  full  =  1  &  Y„  =  0  (40) 

if  P°fty//  =  1  &  Yne  (0,  1] 


for y  =  1 ,...,  8  and  n  =  0,...,  N-  1. 


We  now  comment  on  our  choice  of  solution  ^(8,  p)  for  p  =  0  and  p=  1 .  Our  choice  of 
solution  #8,  p)  for  p  =  0  was  guided  by  the  mathematical  considerations  of  preserving  (a) 
the  continuity  of  g  in  p  at  p  =  0  for  each  8  e  [0,  1  ]  as  well  as  ensuring  (b)  the  continuity  of 
#8,  0)  in  [0,  1],  and  was  reinforced  by  the  circumstance  that  this  choice  causes  no 
conceptual  difficulties:  it  maintains  Ppartiai (y,  pfull)  =  0  for  all  Y  <e  [0,  1]  when  pfull=  0.  Indeed, 
there  are  other  choices  for  the  solution  #5,  0)  which  also  cause  no  conceptual  difficulties 
in  the  above  sense  and  which  maintain  the  continuity  in  (b)  above  but  which  destroy  the 
continuity  in  (a)  above.  An  example  of  this  (and  there  are  many  others)  follows  by  setting 
g(5,  0)  =  74(8)  where,  for  any  fixed  a  e  [0,  1/2),  P£a  is  the  continuous  function  in  [0,  1] 
given  by 
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0 


if  8  e  [0,  a) 


*C(8)  = 


(8  -  a)/(1  -  2a)  if  8  e  [a,1  -  a] 


(41) 


1  if  8  e  (1  -  a,  1]. 

Of  course  the  solution  to  Eqs.  (33)  +  (34)  given  by  modifying  the  second  line  of  Eq.  (37) 
to  g{8,  0)  =  Jfa{8)  lacks  the  continuity  of  g  in  p  at  p  =  0  for  all  8  e  [0,  Vi)  u  (V2,  1].  In  any 

case,  our  actual  choice  of  g(8,  p)  for  p  =  0  is  the  simplest  one,  hence  (for  us)  the  most 
natural. 


Our  choice  of  solution  g(8,  p)  for  p  =  1  was  also  guided  by  the  consideration  of  preserving 
the  continuity  of  g  in  p  at  p  =  1  for  each  8  e  [0,  1  ],  but  this  choice  necessarily  leaves  £(8, 
1)  discontinuous  at  8  =  0  (see  Eq.  (38)).  More  seriously,  this  choice  for  g(8,  1)  implies 
that  Ppartiai(y,  pfuli )  =  pfull  for  ally  e  (0,  1]  when  pfull=  1 ;  but  it  is  not  at  all  clear  that  this  state- 
of-affairs  is  not  physically  valid,  as  follows.  We  observe  that  governing  Eq.  (33),  which  is 
based  upon  the  eminently  reasonable  Eq.  (31a),  becomes,  when  p  =  1  (=  pfull), 

S(Y.1)+S(1  -TU)-S(Y.1)S(1  -  Y,  1 )  =  1  -  (42) 

If  we  now  set  x  =  g(y,  1 )  and  y=  g(  1  -  Y,  1 )  and  consider  the  function  f(x,  y)  =  x  +  y  -  xy 

for  (x,  y)  e  [0,  I]2,  it  is  straightforward  to  see  that  0  <  f(x,  y)  <  1  and  that  f(x,  y)  =  1  iff  (x, 
y)  €  ({1}  x  [0,1])  u  ([0,1]  x  {1}).  In  other  words,  to  satisfy  Eq.  (42)  we  must  have,  for  any 
Y,  either  g(y,  1)  =  1  &  £(1  -  Y,  1)  e  [0,  1],  or  £(1  -  Y,  1)  =  1  and  g{y,  1)  e  [0,  1];  i.e., 
either  fyartial( Y,  1)  =  g{y,  1)  •  1  =  1  &  fyartia,{  1  -Y.  1)  =  g{\  -  Y,  1)  •  1  e  [0,  1],  or  fyartial(  1  - 
Y,  1)  =  1  &  Ppgrtiai  (Y,  1)  e  [0,  1  ]— there  is  no  way  to  satisfy  Eq.  (42)  with  both  fyartia,{Y,  1) 
and  Ppartial{  1  -  Y,  1)  being  less  than  1.  So,  if  we  accept  that  g( y,  1)  is  non-decreasing  in 
Y — as  assumed  following  Eq.  (24) — then  we  must  have  g( y,  1)  =  1  for  Y  e  [V2,  1], 
Alternatively,  we  observe  from  Eq.  (34)  that  we  must  have  ^(0,  1)  =  0  and 
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£(1,  1)  =  1,  and  from  Eq.  (36)  that  we  must  have  g(V 2,  1)  =  1,  so  the  requirement  that 
g{  Y,  p)  be  nondecreasing  in  Y  for  all  p  e[  0,1]  then  dictates  that  g{y,  1)  =  1  for  Y  e  [Vz,  1]. 
So  at  least  for  Y  e  [Vz,  1]  it  must  hold  that  fyartial{y,  pfull)  =  g{ Y,  pull)  pfull=  g{ Y,  1).1  =  1. 


We  have  just  seen  that  g{ 6,  1)  =  1  for  8  e  [Vz,  1],  We  next  consider  8  e  (0,  Vz)-  For  such 
8  we  have  chosen  g(8,  1)  =  1  as  well,  which  choice  implies  the  above  behavior  as  well 
for  Y  e  (0,  Vz),  namely  fyartiai( y,  pfull)  =  1.  This  behavior  is  perhaps  now — in  light  of  the 


previous  argument  for  Y  e  [Vz,  1] — at  least  plausible.  It  is  nevertheless  worthwhile  to 
explore  the  alternatives  to  our  choice  of  g{8, 1 ).  We  categorize  the  possibilities  as  follows: 
either  g(»  ,1 )  is  not  continuous  in(0,  1  ]  or  g{»  ,1 )  is  continuous  in  (0,  1  ].  In  the  first  case, 
one  such  family  of  solutions  for  g(8 ,  1)  is  given  by  setting  g(8,  1)  =  J7a(S)  where,  for  any 
fixed  a  e  [0,  Vz),  Ja  is  the  function  given  by 


.7a  (8)  =  1 


0 

1 


if  8  e  [0,  a) 
if  8  e  [a,  1] 


(43) 


and  therefore 


0 


P'\ 


1 


if  Y  e  [0,  a) 


if  Y  e  [a,  1]. 


(44) 


This  solution  p1^^  (Y,  1 )  exhibits  a  threshold  behavior :  If  Y  <  a  then  pl<Jpartiai (Y,  1 )  =  0  while 

if  Y  >  a  then  p7«art/a/( Y,  1)  =  pull.  The  validity  or  non-validity  of  this  behavior — and 
consequent  appropriateness  of  this  solution — may  easily  be  determined  experimentally. 
In  the  second  case,  in  which  g(»  ,1 )  is  continuous  in  (0,  1  ],  it  turns  out  that  in  fact  the  only 
solution  for  £(-,1)  is  the  one  we  have  already  chosen  in  Eq.  (37).  We  now  demonstrate 
this  claim;  namely: 


Let  g{8,  1 )  satisfy  Eq.  (42)  with  £(0  ,1 )  =  0  and  £(1  ,1 )  =  1  and  let  g{»  ,1 )  be 
non-decreasing  in  [0,  1]  and  continuous  in  (0,  1];  further,  let  £(8,  1)  =  1  for  (*) 
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8  e  [V2,  1],  Then  g(8,  1)  =  1  for  8  e  (0,1]. 

To  prove  this  result  (*)  we  first  make  explicit  a  natural  extension  of  Eq.  (29)  from  basic 
time  segments  to  any  subinterval  of  a  basic  time  segment,  to  wit: 

If  B  =  [mCj1/ 2,  (n+1)xCj1/2)  is  any  basic  time  segment  having  associated 
probability  pfull  =  #U//(B)  and  if  intervals  J'  and  Jare  such  that  J'  c  Jc  B, 

with  Y'  =  Length  {J'}  he, M2  and  Y  =  Length{J}hc,  V2  and  pull[J]  =  Pfun(B)), 

then  we  take 

hariialiy'’  Pfu,l(B))  =  partial  O'' ^  PfJl\)  (45) 

since  the  argument  Y'  inside  the  /^art/a/on  the  LHS  of  Eq.  (45)  refers  to  a 
sub-interval  of  B  of  length  y'tcJ/2  and  the  argument  Y'/Y  inside  the  /^arf/a, 
on  the  RHS  of  Eq.  (45)  refers  to  a  subinterval  of  B  of  length  (y'/y)»Ytc  1/2  = 
y'xCj1/2  as  well.  Further,  we  take 

partial  PfJ M)  =  SW,  PM[Y])pM[Y].  (46) 

The  proof  of  (*)  is  now  as  follows.  Let  A  =  {  8  e  (0,  V2]  \  g(8,  1)  =  1  };  then  V2  e  A  and, 

(t) 

denoting  8*  =  inf  A,  then  8*  e  [0, 1/2].  Now  (8*,  V2]  c  A  :  If  not  then 

(3S#  e  (8*,  V2]){g{8#,  1)  <  1)(tt)  so  that  (V8  €  A)(S#  <  8)(ttt)  [for  if  not  (fff)  then 

(38  e  A)(S#  ^  8)  so  that  g(8#,  1 )  >  g(8 ,  1 )  =  1 ,  a  contradiction  of  (ff)];  hence,  by  (fff),  S# 

is  a  lower  bound  for  A  and  therefore  8*  >  S#  >  0,  a  contradiction  of  S#  >  8*  in  (ff).  Now  if 

8*  =  0  then  we  are  done  since  in  that  case  we  have  by  (f)  that  A  c  (0,  V2],  hence  A  =  (0, 
1/2].  So  suppose  that  8*  >  0.  Since  £(*,1)  is  continuous  in  (0,  1]  then  g(8*,  1)  =  lim5^(5„)+ 

g{8,  1)  =  1 ,  since  (V8  >  8*)(g(8,  1)  =  1  by  (f);  so  fyartial(8*,  pfull)=  g{8*,  1  )pM  = 

1 .1  =  1 .  Denote  pull[8*]  =  fyartial{8*,  pfull)  so  that  pfull[8 *]  =  1 .  Now  S*/2  >  0  and 

g(8*/2,  1 )  =  g( 5*/2,  1 ) .  1  =  £(8*/2,  pJpM  =  Ppartjai(8*l2,  pj  =  ^art/a/((5 */2)/5*,  pM[S*]) 

(47) 
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=  /v/a/(i/2,  Ru„[ 5*])  =  g(te,  null[8*])pfull[8*]  =  g(V2, 1) .  1  =  g(y2,  i)  =  i 

since  Vi  e  A,  where  the  fourth  equality  in  Eq.  (47)  follows  from  Eq.  (45)  with  Y  =  8*/2  and 
y  =  8*  while  the  fifth  equality  in  Eq.  (47)  follows  from  Eq.  (46).  Now  from  Eq.  (47)  we  see 
that  8*/2  e  A;  but  since  872  <  8*  then  8*  *  inf  A,  a  contradiction.  So  it  is  not  true  that 
8*  >  0  and  we  are  done> 


Finally,  another  more  global  example  of  an  alternative  solution  to  Eqs.  (33)  +  (34),  along 
the  lines  of  the  discontinuous  solution  g(8,  1)  presented  in  Eq.  (43),  is: 


r 


g(s,  p)  =  Y 


0 

dvYP) 

i 


V 


if  8  e  [0,  Vi) 

if  8  =  Vz  (pe[0,  1]). 
if  8  e  (V2,  1] 


(48) 


Using  g(8,  p)  we  have,  for  all  pfull  e  [0,  1], 


P*Pa«ia,^Pj=  < 


0 

Pful)  Pfull 


V  P 


full 


if  Y  e  [0,  y2) 

if  Y  =  Vi 
if  Y  e  (V2,  1], 


(49) 


The  question  arises  as  to  where  the  basic  p"  fun[VBW  co]'s  are  to  come  from:  At  this  point 

they  must  be  obtained  externally  to  our  model.  They  may  be  obtained  experimentally  or 
predicted  via  additional  modeling  which  has  to  date  not  yet  been  accomplished.  However, 
once  these  p"  fun[VB M,  co]'s  are  available,  then  the  p" part/a/[ l/EM,  co]'s  may  be  computed  via 
Eq.  (40). 

IV.2.  Individual  y2-Period-Seqmented  SPT  Elements:  Indirect  Probability  Assignment 

As  mentioned  previously  in  Section  111.1 ,  we  feel  that  the  primitive  SPT  elements  are  the 
natural  entities  to  which  disruptive  deviation  probabilities  should  be  assigned;  disruptive 
deviation  probabilities  for  the  y2-period-segmented  SPT  elements  may  then  be  computed 
from  these.  In  this  section  we  present  this  approach.  An  additional  benefit  of  this  approach 
is  that  it  allows  us  to  overcome  the  objection  (alluded  to  parenthetically  following  Eq.  (25)) 
that  the  formulation  of  Eqs.  (24)  and  (28)  for  Ppgrtial  (and  its  elaboration  in  Eq.  (39))  fails  to 
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take  into  account  that  the  Vfe-period-segmented  SPT  elements  sf  have — except  for  j  =1 , 
5 — internal  structure  (Eqs.  (2)  -  (9))  and  therefore 

Pj1  ^31  should  be  based  not  merely  upon  a  single  Yn  for  the  entire  time  segment  but  rather 
upon  as  many  Yn's  for  the  segment  as  there  are  segment  features.  For  example,  s4  should 
use  three  Yn's — one  each  for  the  rising,  constant,  and  falling  time  sub-segments — to 
compute  its  p"  partial  in  accordance  with  whether  the  partial  overlap  of  7rEM[T0,  Tf,  l£M,  co] 
with  s4  occurs  during  the  rising  portion,  the  constant  portion,  or  the  falling  portion  of  that 

segment  or,  more  generally,  during  some  combination  of  two  or  three  of  these.  We 
present  these  enhanced  results  as  well. 

The  reader  may,  upon  first  reading  of  this  paper,  skip  the  details  in  the  present  section 
and  safely  move  on  to  Section  IV. 3;  for  all  that  is  necessary  to  proceed  thereto  is  that 
each  V^-period-segmented  SPT  element  have  a  ;rEM-dependent  disruptive  deviation 
probability  assigned  to  it  and  that  has  been  already  been  accomplished  via  Eqs.  (22)  and 
(40). 

IV.2.1.  Vf  Period-Segmented  SPT  Elements  from  Primitive  SPT  Elements 

In  this  section,  we  derive  expressions  for  the  V^-period-segmented  SPT  elements  in  terms 
of  the  primitive  SPT  elements.  It  is  sufficient  to  consider  in  detail  only  the  zeroth  !/2-period- 
segmented  SPT  elements  sf ;  results  for  sj7,  n  =  1 ,...,  N-  1 ,  then  follow  immediately  from 
Eqs.  (11)  and  (22). 

Consider  firstly  the  four  basic  primitive  SPT  element  types  (A,  B  C,  D),  which  we  label  by 
av,  v  =  A,  B,  C,  D,  and  which  are  given  specifically  by 


aA(f)  =  Mr,  high 

t  e  [0,  xc  1/2) 

(50) 

CTb(0  =  Mr,  low 

t  e  [0,  tCj1/2) 

(51) 

ac(0  =  (A  Mr/ axc)  t  +  Mr,  mid 

t  e  [-axc>1/2,  a xc  1/2) 

(52) 

(7d  (0  =  —  (AMr/pXc)f  +  Mr, 

mid  t  e  [— (3xc  1/2,  pXc1/2). 

(53) 

34 

Approved  for  public  release:  distribution  unlimited. 


A  glance  at  Fig.  2  depicting  the  sf's  reveals  that  each  s°  consists  of  portions  of — but  not 
the  entirety  of  (except  in  the  cases  of  s°  and  s°) — several  c^'s.  These  portions,  which  we 
refer  to  as  partial  primitive  SPT  elements,  are  as  follows: 


aA+  =  aAKaxc,1/2>  xc,1/2) 

(54) 

<*A-  =  cyAf[0,  (1  -  P)xc,i/2) 

(55) 

°A±  =  ^Af[aXC,1/2.  (1  -p)xc,1/2) 

(56) 

aB+  =  f  [Pxc,1/2>  xc,1/2) 

(57) 

aB-=  aBf[0,  (1  -  a)xc  1/2) 

(58) 

aB±=  c?B  I'[Pxc,i/2j  (1  -a)xc,i/2) 

(59) 

^C+^cK0’  «xc,1/2) 

(60) 

^C“=CfcK-axc,1/2,  0) 

(61) 

aD+  =  aD  ftO,  Pxc,1/2) 

(62) 

^D-=  aD  l'[-Pxc,1/2i  0); 

(63) 

further,  it  is  not  ac-and  aD-that  are  directly  needed  but  rather  their  time  translates 


°C-(0  =  ^C-(^-xc,1/2) 

te  [(1  -  a)xc,1/2>  xc,1/2) 

(64) 

°D-(0  =  CTd-  (f  —  Xc,i/2) 

t  e  [(1  -  P)xC,1/2>  xc,1/2)- 

(65) 

We  denote  by  aM,  p  e  {A+,  A-,  A+,  B+,  B-,  B±,  C+,  C-,  D+,  D-}  =  M,  any  of  the  various 
partial  primitive  SPT  elements  indicated  in  Eqs.  (54)  -  (65)  (but  not  in  Eqs.  (61)  and  (63)) 
and,  further,  by  1^  the  respective  restriction  interval  associated  with  cr^  (e.g.,  Ic-  =  [(1  - 
oc)xCji/2j  tc,i/2))-  We  have  written  the  c^'s  in  the  above  as  restrictions  to  subintervals  of  the 
full  domains  of  definition  of  the  basic  primitive  SPT  elements  in  order  to  emphasize  that 
the  disruptive  deviation  probabilities  to  be  assigned,  in  the  next  subsection,  to  the  aM's  for 
H  e  9A.  are  in  fact  /^art/a,( y,  pfu//)'s  in  the  sense  of  the  previous  subsection,  where  pfull  is  that 
of  the  “parent”  basic  primitive  SPT  element  av(|a)  of  (e.g.,  aA  is  the  parent  of  each  of 
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aA+,  aA-,  and  aA±  so  that  if  jlx  =  A+,  A-  or  A±  then  v(n)  =  A)  and  Y  is  the  relative  fractional 
length  given  by 

Y  =  %=  Length^} /Length{“parent”  ov((l)}.  (66) 

In  this  connection,  we  also  observe  that  if  we  denote  <fc (t)  =  oc(t-xc>V2) for  t  e 

[(1  -  a)xc  1/2,  (1  +  cc)xCji/2),  and  <j*D(t)  =  oD(t-  xc  1/2)  for  t  e  [(1  -  p)xci/2,  (1  +  P)xc  1/2),  then  in 
terms  of  these  time  translates  we  have 

<*C~  =  C^cKC1  -  a)Tc,1/2>  *0,1/2)  and  °D-  =  C5of[(1  -  P He, 1/2)  *<5,1/2)-  (67) 

This  observation  will  be  relevant  in  assigning,  in  the  next  subsection,  disruptive  deviation 
probabilities  to  ac-  and  cjd-,  for  we  will  continue  to  assume  for  primitive  SPT  elements 
(both  full  and  partial) — as  we  have  previously  assumed  for  Vb-period-segmented  SPT 
elements  (see  Eq.  (22)) — the  property  of  time  translation  invariance  (or  what  we  have 
previously  termed  “modularity”)  of  disruption  probabilities. 

We  denote  by  a°  the  extension  of  to  [0,  xc>1/2),  i.e., 

r  c^(0  if  f  G  b 

aj(f)  =  ^  (n  e  94.)  (68) 

lo  if  f  e  [0,  Te)  \  I,. 

Since  aA  and  aB  in  their  entirety  are  required  to  construct  1/2-period-segmented  SPT 
elements  s°  and  s5°,  we  also  allow  the  index  p  to  take  on  values  A  and  B  (but  not  C  or 
D),  with  /A  =  [0,  xc  1/2)  =  Iq ;  and  we  write  94+  =  94  u  {A,  B}.  Thus  the  specification 
jli  g  in  Eq.  (68)  is  to  be  replaced  by  n  e  9\i+.  Note  that  aA  =  aA  and  ctb  =  aB. 


At  this  stage,  all  V^-period-segmented  SPT  elements  in  [0,  xc  1/2)  may  be  constructed  as 
sums  of  the  (full  or  partial)  primitive  SPT  elements  of  Eq.  (68)  (with  p  e  94+).  To  wit,  for 
all  t  e  [0,  xCii/2)  we  have: 


^0  o 
Si  =  aA 


Oo  o  o 
S2  -  ac+  +  aA. 


o0  0  0 
S3  _  °A  +  aD 
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(70) 


„o  o 

s5  -  Ob 


„o  0  0 

S6  _  aD+  +  °B 


„0  0  0 
S7  -  °B-  +  0C- 


sS 


0,0,0 
aD+  +  °B±  +  ac-  ■ 


We  denote,  for y  =  1,  8, 

=  {(a  e  fM+ 1  a°  occurs  in  the  sum  for  s°} 
so,  for  example,  =  {  A+,  C+,  D-};  also 

(71) 

IV.2.2.  Primitive  SPT  Element  Disruptive  Deviation  Probabilities 

We  next  compute  disruptive  deviation  probabilities  for  the  partial  primitive  SPT  elements 
based  upon  disruptive  deviation  probabilities  for  the  basic  primitive  SPT  elements. 

To  begin,  we  denote  by  pv  full[ \/EM,  co],  v  =  A,  B,  C,  D,  the  probability  that  basic  SPT  element 
cv  suffers  a  disruptive  deviation  when  it,  or  a  time  translate  of  it,  occurs  as  part  of  an  N- 
element  SPT  and  when,  in  addition,  it  is  fully  encompassed  by  some 
xEM[T0,  Tf,  \/EM,  co].  As  previously  assumed  for  V^-period-segmented  SPT  elements,  we 

again  assume  modularity  for  these  basic  primitive  SPT  elements:  The  probability  of  a 
disruptive  deviation  of  each  one  of  the  four  basic  SPT  primitive  elements  is  independent 
of  where  that  primitive  element  occurs  in  the  /V-element  SPT,  as  long  as  that  primitive 
element  is  completely  encompassed  by  x™.  And  once  again,  the 
Py  fUJVEM’  c°]'s  must  Pe  obtained  externally  and  made  available  to  our  model ;  but  note 

that 

PA,  full  '■'EM’  0)]  =  Pi°  full^EW  0)]  and  PB  fui^EW  0)]  =  P5  fuii^E M>  ro]-  (7^) 
Also,  consistently  with  our  previous  notational  convention,  we  denote  for  convenience  p° 
J  ^EM’  °>]  =  PA,  fJ  VEW  °>]  and  PB°  fuil  ^EM’  ®]  =  PB,  full  VEW  ®]- 


In  addition  to  the  four  probabilities  pv  fujl[VEM,  co]  above,  we  also  need  expressions  for  the 
disruptive  deviation  probabilities,  p°[  VEM,  co],  of  the  g°  ,  |i  e  IM .  Consider  first  y  =  A+:  Since 

the  parent  of  ga+  is  ga  then  from  Eq.  (66)  we  have  L  =  1  -  a,  thus,  using  Eq.  (40), 
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Pa°+[Kem.  ®]  =  Ffaartiai ( 1  ~  «,  p°  M)  =  1  -  (1  -  p°  "  for  p°  M  e  [0,  1]  &  0  6(0,1) 

(73) 

(where  we  have  taken  into  account  that  a  =  0,  1  is  precluded  by  our  initial  definition  [in 
Section  111.1]  of  the  type-C  SPT  primitive  element).  Similarly,  we  have  JA-  =  1  -  p  and 


PA0i'/EM>®]=  ^a/(1-P.PA°M)=1-(1-Pfl°J1_P  f0r  PA.fti//  e  11  &  P^0’1) 

(74) 

(also  where  p  =  0,  1  was  precluded  earlier).  Next  we  have  yA±  =  1  -  (a  +  p)  and 

<o]  =  -(o  +  P).Pft,M)  =  1  - 

forpA0M,s[0,  1]  &  a  +  pe(0,  1)  (75) 

(since  a  +  p  =  0,  1  has  also  been  precluded  previously).  The  expressions  fo rpB°+[\/EM,  to], 
Pb-EKem’  ®],  and  PB°+EK =m>  ®]  follow  directly  from  those  for  p°+[VEM,  co],  pA° _[VEM,  co],  and 
p°+[VEM,  co],  respectively,  via  the  replacement  everywhere  in  the  latter  three  expressions 

of  pA  full  by  pB  full  as  well  as  the  interchange  everywhere  of  a  and  p  since  yB+  =  1  -  p,  v 

=  1  -  a,  and  yB±  =  i  -  (a  +  P).  Finally,  since  yc+  =  Vz  =  yc~,  with  the  parents  of  ac+  and  ac- 

being  ac  or  time  translate  <j*c  respectively;  and  yD+  =  Vz  =  yD~,  with  parents  of  aD+  and  gd- 
being  aD  or  a*D  respectively;  then  we  have  from  Eq.  (40)  that 


Pc+E^EM’  ®]  _  Pc-[^EM’  ro]  -  1  ("I  Pc,  full )  ^  Pc,  full  e  [0> 

pSSVe m>  ®]  =  Pq-E^em’  co]  =  1  -  (1  -  pD  J1/2  if  PD  full  e  [0,  1], 


(76) 


The  computation  of  the  disruptive  deviation  probabilities  for  the  partial  primitive  SPT 
elements  is  now  complete. 

IV.2.3.  VrPeriod-Seg merited  SPT  Disruptive  Deviation  Probabilities:  Full  Overlap 
of  nEM  and  sj1 
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We  now  compute  the  disruptive  deviation  probabilities,  pf  M[\/EM,  to],  y  =  1 . 8,  for  the 

1/2-period-segmented  SPT  elements,  s°,  using  the  disruptive  deviation  probabilities 
pf[VEM,  to],  jii  €  9\4+,  of  the  primitive  SPT  elements  cr°,  \x  e  9\i+,  under  the  condition,  as  in 
Section  IV.1 .1 ,  that  there  is  full  overlap  of  nEU  and  sf.  We  suppress  the  argument 
[VEM,  to]  of  pf.  Using  Eq.  (69)  along  with  reasoning  analogous  to  that  presented  directly 
following  Eqs.  (31)  we  find 


P,°  M  =  PA°  fun  (77) 

P2°M=  PC>(1-PC0;PA°.  (78) 

P3°  u  =  PA-  +  0  -  PA°-)  ft°-  (79) 

Pi”  full  =  ftV  (1  -  ft°J  PA°±  +  (1  -  ft°+)(1  -  PA°±)  ft°-  (80) 

P5°  full  =  Pb,  full  <81) 

Pe°  full  =  PD°+  +  (1  -  PD0,)  P°+  (82) 

P7°m  =  Pb0-  +  (1  -Pb°-)P°-  (83) 

Ps°  full  =  PD°+  +  (1  -  PD°+)  PB±  +  (1  -  P0°J(1  -  PB±)PC°-  (84) 


IV.2.4.  VrPeriod-Segmented  SPT  Disruptive  Deviation  Probabilities:  Partial 
Overlap  of  nEU  and  s f 

Finally,  we  compute  the  disruptive  deviation  probabilities,  p°part(a/[l4M.  00],  y  =  1,  8,  for 

the  !/2-period-segmented  SPT  elements  sf,  using  the  disruptive  deviation  probabilities, 
pf[VEM,  co],  |ii  g  9\i+,  of  the  primitive  SPT  elements  <r°  \x  e  M+,  under  the  condition,  as  in 
Section  IV. 1.2,  that  there  is  partial  overlap  of  7rEM[T0,  Tf,  VEM,  co]  and  sf.  We  first  define 
(see  Eq.  (23)),  for  p  g  M, 

V  JJo,  Tf)  =  Length{[T0,  Tf)  n  IJ/ Length^}.  (85) 
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Then — still  using  the  assumption  that  P  ial  depends  only  upon  Y  (however,  see  below) — 


we  have 


Pi  ,  partial  Ppartial  (  V  A’  PA,  fuJ  (86) 

P2%,«  =  Ppa^IVc^PcV  +  l1  -Ppa„J\0,'PcJIPparlJ\A,’P°J  (87) 


P3%a«  =  Ppa«0'„,flJ  P°_)  +  [1  -p^P P^JlP^Po^,  p»_) 


P4°pa««  =  P^Po.Ca.P^  +  11  "  P^Po,  C*.  P^lP^Po,  PA° J 

+  [’-  'WA  c..  Pc°+)1 11  -  Pp.«  Po.  a*  PA°t)]  Pp^Pa  D-.  PD°-) 

P5. partial  =  Ppaniai Po,  B>  PB(  M) 

P6°  pa«  =  PpppppPo.  D+.  f&)  +  [1  -  Ppa„Po,  D+’  PD°+)1Ppa«(T0.  Bt  PB+> 
P/°  pa«  -  Ppa„MPo,  B-.  PB°_)  +  M  -  Ppp„MPo,  B-.  PB°J]Ppa«P„,  c-  Pc°_) 
Ar  partial  =  P^Po.  Dt  PD°.)  +  [1  "  Ppa« P'0.  Dt.  PD°J) Pp,* Po,  B±.  PB°±) 

+  I1  -  Ppa„a)Po,  Df  Pnt>]  I1  -  P„,M  Po.  Bf  PrJI  P„,„jr 0,  0-.  P°J- 


(88) 

(89) 

(90) 

(91) 

(92) 

(93) 


partial 


B±/J  ^ partial 


In  summary,  we  have  computed,  for y  =  1 ,  ...,  8,  both  p°  fa//[\/EM,  ®] — wa  Eqs.  (77)  through 
(84)— and  p°  partial  Ivbw  wa  Eqs-  (86)  through  (93)— from  basic 

primitive  SPT  elements  <\,  v  =  A,  B,  C,  D,  with  the  aid  of  the  intermediary  partial  primitive 
SPT  element  disruptive  deviation  probabilities  P°[VEM,  co],  |u  e  M. 


The  assumption  that  ppartjal{ Y0,  P°,  M)  depends  only  upon  its  two  listed  arguments  may 

be  simplistic  in  principle  in  case  |u  =  C+,  C-,  D+,  D-  since — on  the  one  hand — it  is 
reasonable  to  conjecture  that  the  (notationally  suppressed)  voltage  dependencies  of  the 
various  P  .J s  (as  well  as  of  the  various  Pfull's)  are  functions  of  the  ratio  VEM  /\4ignai  puise 
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and — on  the  other  hand — it  is  the  case  that  the  partial  primitive  SPT  elements  associated 
with  the  above  four  p's  do  not  have  constant-in-time  signal  voltage  levels.  Indeed,  for 
each  of  these  four  partial  primitive  SPT  elements,  the  I4ignai  pU|Se  voltage  range 

encompassed  during  any  sub-intervals  thereof  depends  not  only  upon  the  length  Y0  p(7^, 
Tf)  of  the  sub-interval  but  also  upon  the  position  of  that  sub-interval  within  the  element. 
Since  differing  such  intervals  may  be  completely  characterized,  hence  distinguished,  by 
giving — in  addition  to  their  length — also  the  location  of  the  interval’s  center  point,  say  x0 

n(4  Tf)  =  Midpoint  of{[T0,  Tf)  n  4},  then  x0  „  may  augment  the  argument  list  of  j^art/a/(Y0 
p°  full)  to  yield  P#partia\ Y0  P°  M)  and  thus  provide  an  “anchor”  to  specify  interval-location- 
dependent  quantities,  e.g.,  VEM  / \4ignai  pulse  at  To  ^  uPon  which  Ppartial  may  be  declared  to 
depend  in  some  model;  of  course  this  dependence  must  further  be  specified.  Eqs.  (87)  - 
(89)  and  (91)  -  (93)  may  then  be  modified  by  replacing,  on  their  RHS's,/^art/a/(Y0  p,  p°  full), 

p  =  C+,  C-,  D+,  D-,  with  Ppartiat (Y0  p°  full) .  We  will  not  pursue  this  refinement  further  in 
this  paper  and  take  Eqs.  (86)  -  (93)  as  they  stand  as  our  model  for  p°  partjal [ \/EM,  co],  j=  1 , 

...,8. 

IV.3.  Arbitrary  A/-element  SPT 

In  the  previous  sections  we  have  shown  how  to  build  an  A/-segment  SPT  from  Vfe-period- 
segmented  SPT  elements  and,  further,  how  to  assign  (in  two  ways)  disruptive  deviation 
probabilities  to  these  elements.  We  now  combine  these  two  to  compute  the  disruptive 
deviation  probability  for  the  entire  idealized,  normal,  /V-segment  SPT  To  this  end,  and 

with  notation  as  in  Section  III,  fix  7rEM[T0,  Tf,  \/EM,  co]  and  consider  any 

e  n7fN[0,  Te],  with  jn  e  {1 . 8}  for  n=  0,...,  N-  1,  and  write 

n71  =  E^To  s fn  with  (Vn  =  0,  ...,  N-  1  )(<  jn,  jn+ 1  >  e  J* )  as  in  Eq.  (14).  Now  consider  the 

following  conceptual  random  experiment:  At  t  =  0,  irradiate  with  nEM  a  single  signal  line 

carrying  n71  and  observe  the  first  time  segment  in  which  a  disruptive  deviation  of  the 

SPT — via  disruptive  deviation  of  its  l/2-period-segmented  SPT  element  in  that  time 
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segment — takes  place  (which  disruptive  deviation  can  be  identified  by  the  mal-execution 
of  the  jl.iC  task,  that  is  executing  and  to  which  a*  belongs).  In  this  situation,  it  is  not  all 

of  n71  that  is  of  interest  but  rather  n^Wn71  n  ^EM[.  The  possible  outcomes  of  the  random 
experiment  can  be  given  by  the  members  of  the  set 

Tflnmin+J[—nmax-^nmJT0’  Tl  "max*  1  )>  (94) 


with  the  outcome  “n”  representing  the  situation,  when  n  <  nmax,  that  the  first  time  segment 
in  which  a  disruptive  deviation  of  the  SPT  takes  place  is  n  and  with  the  outcome  “nmax+ 
1”  representing  the  situation  in  which  no  disruptive  deviation  takes  place  in  any  time 
segment  occurring  during  [7^,  Tf).  The  set  of  events  associated  with  QPb’  V  is  of  course 

n[T0,  Tf) 

taken  to  be  2  N  (the  power  set  of  QF6’  V).  We  now  specify  a  probability  measure, 

o[t0,  V 

denoted  <PjJ o-  Tt\  on  (Q^o-  V,  2  N  ).  First  we  take,  for  n  =  nmjn,...,  nmax, 


Vp  r'>(M)  =  Pjj  VEU,  to;  To,  7;]  -  P\  [  VEU,  to;  T„,  7J])  (95) 


where 


Pj,  full[  ^EM’  ®] 


ifr=n™n+1--->nmax-1 


Pj [  T0,  Tf] 


(96) 


Pj,  partial  ^EM’  031  ^0>  Tf]  if  -  ^mjn  Or  H 


max 


and  possibly  p|,  part/a,[  VEM,  co;  T0,  Tf]  =  fu//[\/EM,  co]  when  r  =  nmjn  or  nmax;  clearly, 


0  <  (Pjjo’  Tf]{{n})  <  1 .  Note  that  p-,  partiai  may  be  that  from  either  of  the  two  alternative 
formulations  given  in  Section  IV. 2.  It  is  straightforward  to  show  that 

r'l«n»  = 1  -  n,1™i  (1  -  p?j  *  1  (97) 

and  we  therefore  have 

1 » =  n  (i  -  fy,  os) 
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as  expected.  Since  {{nH  n_"lax+1  is  a  partition  of  Q[7b’  V  we  may  then  extend  <P[7o -  V  to  all 

“  n-  nmin  r  N  3  N 

nlT0,  Tf )  [T0,  Tf) 

of  2  N  by  0’  Tf\A)  =  Y.neA  Pjjo’  ^({n})  for  A  e  2  N  ;  then  Pjjo,  Tf),  being  manifestly 

countably  additive,  is  thus  a  probability  measure.  Our  probability  space  for  observing  the 
occurrence  of  a  disruptive  deviation  induced  by  7tem[T0,  Tf,  VEM,  co]  upon  a  single  signal 

line  carrying  a  specific  but  arbitrary 

n[T0,  Tf) 

77^[0,  Te]\]7t^  n  ^-EM[  is  then  (Q jjo-  Tf\  2  N  ,  PjJ0,  rf)^  where  we  have  denoted  by 
/7*[0,  Te]  \]i?  n  tiem[  the  set 

n\[Q,  te]\]^  n  ^EM[  =  {^*f  n  ;rEM[  |**e77*[0,  TE]}.  (99) 

The  event  “SPT  n71  suffers  a  disruptive  deviation”  is  <EV  D  =  V  \  {  nmax+  1}  with, 
according  to  Eq.  (98), 

V (<£,.  o)  =  1  -  n  q"™L  d  -  p?q  [  ''em.  To,  7J])S1.  (100) 

This  is  our  desired  final  expression  for  the  probability  <P^>-  ^)(^d)  that  SPT  ^suffers  a 
disruptive  deviation  when  subject  to  xEU. 

It  is  interesting  to  observe  that  the  result  given  in  Eq.  (100)  can  be  obtained  in  an 
alternative  fashion  using  an  ( nmax  -  nmjn  +  1)-factor  product  space  formulation  in  which  a 

factor  in  the  product  appears  for  each  of  the  ( nmax  -  nmjn  +  1)  V^-period-segmented 
elements  of  the  SPT  restricted  to  \7t>t  n  ;rEM[  and  in  which  the  two  possible  outcomes 

achievable  for  each  individual  y2-period-segmented  element — a  disruptive  deviation  of 
that  element  occurs  or  no  disruptive  deviation  of  that  element  occurs — are  taken  to  be 
independent  between  elements.  That  is,  if  for  n  =  rimjn, nmax\Ne  take  Gn  =  {d,  ~c/},  where 

d  is  the  outcome  for  s £  that  s fn  suffers  a  disruptive  deviation,  with  fn{{d})  =  p f  ,  and 

fn{{~d})  =  1  -  p",  and  set 


43 


Approved  for  public  release:  distribution  unlimited. 


and 


©'r°’  V  =  Xn"™x 

N  ''-'win 


(101) 


V(x"T-  A, )  =  x 

A/  1 1  —  1 1  mi  n 


max 
n  =  nmir 


(An)  for  Ane?\ 


(102) 


then  in  this  setting  the  event  “SPT  **  does  not  suffer  a  disruptive  deviation”  is  given  by 


«[J0’  Tf) 


~(Z>,  D) 


=  {<~d,  ~d,  ...,  ~cb>}  e  2Pw  with  probability  fjjo-  D))  =  Tlqn™*jn{1  -  pty, 


which  probability  is  consistent  with  the  result  in  Eq.  (100).  More  generally,  the  event  {/?} 
in  the  first  formulation  corresponds  to  the  event  of  the  second  formulation  given  by 

Sfo  D)  n=  {<~cf, ~  d,  d,  o,...,  o>  e  ©Ob’  V  |  ~d occurs  in  the  first  n-  1  entries 


&  o  e  {d,  ~d}} 

=  h)xmx(x*>4,  e„);  d03) 

hence, 

where  the  final  equality  follows  from  Eq.  (95). 


The  reason  we  have  not  used  this  second,  product  space  formulation  is  because  many 
of  the  outcomes  in  ©Ob’  Tf]  are  not  appropriate  for  our  physical  context.  For  example,  an 

outcome  in  ©0b’  7?)  of  the  type<cf,  ~d,  ...,~d> — in  which  a  “no  disruptive  deviation”  outcome 

must  be  observed  for  every  individual  l/2-period-segmented  SPT  element  subsequent  to 
the  first,  with  each  such  subsequent  element  following  an  individual  element  (the  first 
element)  for  which  a  “yes  disruptive  deviation”  outcome  has  occurred — may  not  be  an 
observable  outcome  at  all  since  the  SPT  may  cease  to  exist  after  its  first  element  is 
affected;  further,  the  product  formulation  in  general  ascribes  a  positive  probability  to  such 
outcomes  (similarly  to  Eq.  (104)).  In  order  to  avoid  considering  such  perhaps 
unobservable  outcomes,  we  have  chosen  our  outcomes  to  be  of  the  type  given  in  our  first 
formulation,  namely  {n},  in  which  observation  of  elements  temporally  beyond  the  first 
affected  one  is  not  required. 
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V.  A  Random  Collection  of  SPT’s 

In  the  above,  we  formulated  a  probabilistic  model  for  the  disruptive  deviation  of  an 
arbitrary  but  fixed  e  77^[0,  Te]  that  is  being  carried  on  a  given  signal  line  when  that 

line  is  exposed  to  EM  disturbance  xEM[T0,  Tf,  VEM,  co],  However,  as  discussed  in  the 
Introduction  to  this  paper,  in  order  to  determine  the  probability  of  upset  by  ^EM  of  a  |uC 
executing  a  particular  application  task  a  e  ?4,  we  must  resort  to  analyzing  a  stochastic 

mix  of  piC  tasks  %  =  <[\a]#,  [df,  e  C*  x  i#xD i.e.,  to  analyzing — for  each  of 

the  piC  signal  lines — a  stochastic  mix  of  n71'^  on  that  signal  line.  Since  in  this  paper  we 
have  been  considering  only  a  single  such  signal  line,  then  we  must  therefore  consider 
analyzing  a  stochastic  mix  of  7^'s  on  that  signal  line.  Thus  we  next  want  to  allow  the 

possibility  that  the  specific  SPT  on  that  line  during  xEU — that  SPT  being  any  one  of  the 
possibilities  k%\\k%  n  /z-EM[  e  77^[0,  TE]f]^  n  ^EM[ — is  determined  stochastically  (for 

example,  each  of  the  possibilities  for  ^EM[  on  the  line  is  equally  likely.)  We  may 

easily  write  the  expression  for  the  disruptive  deviation  probability  in  this  case  as  a 
superposition  of  the  SPT  results  in  the  previous  Section  IV. 3. 

To  that  end,  we  first  abbreviate  the  notations  for  T^Wn71  n  7tem[  and 

77^[0,  Te]  Wn71  n  7tem[  to  n71  f  and  77^[0,  TE]  \  respectively.  Next,  we  index  the  members 

jP\ of  77*  [0,  Te]\,  writing  ^sf,  s=  1,  ...,  #(77*  [0,  TE]  [)  =  S for  them,  where 

#(77*  [0,  Te]\)  =  2nmax~nmin  +  3.  We  then  assign  probabilities  P  W/r*’s[),  s=  1,  ...,  S, 

"  A/I 

such  that 

Zss.,4*|<^si>  =  1-  (105> 

This  assignment  reflects  our  knowledge  of  the  probabilities  of  occurrence  of  the  various 
possible  /r*[’s  on  the  signal  line,  which  knowledge  results  from  a  random  experiment 
conducted  on  the  line  in  the  absence  of  any  EM  disturbances  and  for  which  the  outcome 
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set  is  n  ^[0,  Te]  f  or,  alternatively,  from  some  model  of  the  relative  frequencies  of 
occurrence  of  all  possible  j^\ s  on  the  signal  line.  For  example,  equal  probabilities  for  all 
of  the  n71  f’s  yields  Pnv^x*’s\)  =  1/S,  s  =  1 ,  S.  Next,  the  outcome  set  for  a  different 

random  experiment — an  experiment  to  determine  the  probability  of  disruptive  deviation 
on  the  stochastically  populated  line — is  simply  Q;/^  =  {D,  ~D},  where  “D”  is  the  outcome 

that  the  signal  line  suffers  a  disruptive  deviation  (as  opposed  to  the  previously  used 
outcome  “d”  to  signify  that  a  single  V^-period-segmented  SPT  element  suffers  a  disruptive 
deviation).  Finally,  we  assign  a  probability,^ ’  W{D}),  to  {D}.  Using  Eq.  (100) 

to  express  Tf)(‘EVD)  for  any  A/-element  SPT,  with  this  symbol  augmented  in  notation 
by  a  subscript  “s”  to  indicate  any  particular  SPT  /r^sf  to  which  it  refers  (i.e.,  we  write 
(P  i  ft  7/)('£p,D)),  we  then  take 

P„*j({D})  =  If,,  s  1  (106) 

(and  of  course  ^^({-D})  =  1 

This  is  our  desired  final  expression  for  the  probability  P  W{D})  of  disruptive  deviation  of 

A/I 

an  SPT  on  the  signal  line  in  the  stochastic  mix  case. 

VI.  Conclusion 

We  have  presented  a  model  that  provides  predictions  of  the  occurrence  of  disruptive 
deviations  of  signal  line  activity  in  a  jliC  executing  an  application  task  when  the  |uC  is 
subjected  to  external  irradiation  by  an  EM  pulse  impinging  upon  it,  modeling  this  situation 
from  a  signal-centric  point  of  view.  We  have  specified  the  set  of  all  possible  idealized, 
normal,  A/-element  SPT’s  that  may  be  carried  on  any  signal  line  of  the  piC  and  have 
described  as  well  the  interaction  of  any  such  SPT  with  a  rectangular-envelope-modulated 
monochromatic  sinusoidal  pulse  influencing  the  SPT.  This  interaction  may  result  in  a 
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disruptive  deviation  of  any  SPT  on  the  line  present  during  the  EM  pulse  time  window  and, 
consequently,  in  disruption  of  jliC  task  execution — an  “upset”.  Upsets  occur  stochastically 
and  we  have  provided  expressions  for  the  probability  of  such  upsets,  both  for  any 
specified  possible  SPT  on  a  single  signal  line  as  well  as  for  stochastic  collections  of  such 
SPT’s  on  that  single  line.  (Results  for  collections  of  such  signal  lines  [with  a  complete  piC 
being  a  very  large  such  collection]  will  be  presented  in  a  subsequent  paper.)  The 
disruptive  deviation  probability  for  an  A/-element  SPT  is  given  in  terms  of  disruptive 
deviation  probabilities  for  each  of  its  N  elements;  these  element  disruptive  deviation 
probabilities  are  not  calculated  in  the  model  but  rather  must  be  provided  externally,  either 
by  being  obtained  experimentally  or  via  additional  modeling  performed  beyond  the 
present  model  and  which  has  to  date  not  yet  been  accomplished.  Further,  the  disruption 
probability  associated  with  any  temporal  sub-interval  of  any  V^-period-segmented  SPT 
element,  which  subinterval  is  also  completely  overlapped  by  the  EM  pulse  time  window, 
is  taken  most  simply  to  be  a  function  of  the  length  of  that  sub-interval  as  well  as  of  the 
disruptive  deviation  probability  of  the  complete  V^-period-segmented  SPT  element  in 
which  it  resides.  At  a  more  detailed  level,  the  sub-interval’s  disruption  probability  takes 
into  account  not  only  its  length  but  also  the  particular  structure  of  the  host  V2-period- 
segmented  SPT  element  itself. 


References 

Clarke,  T.,  Taylor,  A.,  Estep,  N.,  Yakura,  S.,  Brumit,  D.,  Dietz,  D.,  Hemmady,  S.,  Duffey, 
J.,  "Predictive  modeling  of  high-power  electromagnetic  effects  on  electronics", 
Electromagnetics  in  Advanced  Applications  (ICEAA),  201 1  International  Conference,  pp. 
429-431,  12-16  Sept.  2011,  doi:  10.1 109/ICEAA. 201 1 .6046377. 

Clarke,  T.,  Dietz,  D..  French,  D.,  "Modeling  of  radio-frequency  effects  on  a 
microcontroller",  Electromagnetics  in  Advanced  Applications  (ICEAA),  2012  International 
Conference,  pp.1266-1269,  2-7  Sept.  2012,  doi:  10.1 109/ICEAA.201 2.6328823. 

Estep,  N.A.;  Petrosky,  J.C.;  McClory,  J.W.;  Kim,  Y.;  Terzuoli,  A.J.,  "Electromagnetic 
Interference  and  Ionizing  Radiation  Effects  on  CMOS  Devices,"  Plasma  Science,  IEEE 

47 


Approved  for  public  release:  distribution  unlimited. 


Transactions  on,  vol.40,  no. 6,  pp.1 495-1 501 ,  June  2012,  doi: 
1 0.1 109/TPS.201 2.21 93600. 

Estep,  N.  A.  Combined  Effects  of  Electromagnetic  Interference  (EMI)  and  Ionizing 
Radiation  on  Digital  Inverters.  Nuclear  and  Space  Radiation  Effects  Conference 
(NSREC),  July  2010. 

Estep,  N.A.  High  Power  Microwave  (HPM)  and  Ionizing  Radiation  Effects  on  CMOS 
Devices.  Master’s  Thesis,  Air  Force  Institute  of  Technology,  Mar  2010. 

Firestone,  T.  M.  RF  Induced  Nonlinear  Effects  in  High-speed  Electronics,  MS  Thesis, 
University  of  Maryland,  2004. 

Firestone,  T.  M.,  J.  Rodgers  and  V.  L.  Granatstein.  Response  of  CMOS  Logic  with 
Electrostatic  Discharge  Protection  to  Pulsed  Microwave  Excitation.  Submitted  to  IEEE 
Transactions  on  Circuits  and  Systems,  2006 

Granatstein,  V.  L.  Effects  of  High  Power  Microwave  and  Chaos  in  21st  Century  Analog 
and  Digital  Electronics.  MURI  Final  Performance  Report,  AFOSR  Grant  Number 
F4962001 10374,  2007. 

Flemmady,  S.,  X.  Zheng,  E.  Ott,  T.  M.  Antonsen  and  S.  M.  Anlage.  Universal  Statistics  of 
the  Scattering  Coefficient  of  Chaotic  Microwave  Cavities.  Phys.  Rev.  E  71 , 056215,  2005. 

Henderson,  R.;  McMasters,  D.;  French,  D.;  Clarke,  T.,  "Direct  power  injection  of 
microcontrollers  in  PCB  environments",  Electromagnetics  in  Advanced  Applications 
(ICEAA),  2012  International  Conference,  pp.1 149-1 151 ,  2-7  Sept.  2012,  doi: 
1 0.1 1 09/ICEAA.201 2.6328803. 

Laurin,  J.,  S.  G.  Zaky  and  K.  G.  Balmain.  On  the  Prediction  of  Digital  Circuit  Susceptibility 
to  Radiated  EMI.  IEEE  Transactions  on  Electromagnetic  Compatibility,  Vol.  37  No  4,  528- 
535,  1995. 

Liu,  C.  K.  and  C.  J.  Ho.  Frequency  Dependence  of  AC-Interference  Susceptibility  in  Digital 
ICs.  IEEE  Proceedings  on  Science  Measurement  &  Technology.  Vol.  142,  2,  183-18, 
1995 

48 


Approved  for  public  release:  distribution  unlimited. 


Richardson,  R.  E.  Modeling  of  Low-Level  Rectification  RFI  in  Bipolar  Circuitry.  IEEE 
Transactions  on  Electromagnetic  Compatibility,  Vol.  21, 307-311,  1979. 


Taylor  A.  Microcontroller  (8051 -Core)  Instruction  Susceptibility  to  Intentional 
Electromagnetic  Interference  (IEMI).  To  be  submitted  as  Master’s  thesis,  University  of 
New  Mexico,  201 1 . 

Vick,  R.  and  E.  Habiger.  The  Dependence  of  the  Immunity  of  Digital  Equipment  on  the 
Hardware  and  Software  Structure.  Proceedings  of  International  Symposium  on 
Electromagnetic  Susceptibility,  ed.  L.  Zhang  and  K.  Zhou,  383-386,  1997. 

Wang,  H.,  S.  V.  Rodriguez,  C.  Dirik  and  B.  Jacob.  Electromagnetic  Interference  and 
Digital  Circuits:  An  Initial  Study  of  Clock  Networks.  Electromagnetics,  26,  73-86,  2006. 

Wendsche,  S.,  R.  Vick  and  E.  Habiger.  Modeling  and  Testing  of  Immunity  of 
Computerized  Equipment  to  Fast  Electrical  Transients.  IEEE  Transactions  on 
Electromagnetic  Compatibility,  Vol.  41  No  4,  452-459,  1999. 

Yang,  H.  Y.  D.  and  R.  Kollman.  Analysis  of  High-Power  RF  Interference  on  Digital  Circuits. 
Electromagnetics,  26,  87-102,  2006. 

Acknowledgement 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research. 


49 


Approved  for  public  release:  distribution  unlimited. 


Amplitude 


tlx 


C,  Vz 


50 


Approved  for  public  release:  distribution  unlimited. 


Amplitude 


Amplitude 


Amplitude 


Amplitude 


Amplitude  Amplitude  Amplitude  Amplitude 


Figure  2.  The  eight  V^-period-segmented  SPT  element  types  (n  =  0). 
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Figure  3.  An  eight-segment  SPT  constructed  from  V^-period-segmented  SPT  elements. 
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